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This  report  is  diviae#*irTwTwo  main  parts:  Experimental  Work  (B)  and 


Theoretical  Work  (C). 

In  the  experimental  part  we  first  report  the  noise  in  gold  metal  films. 

1  2 

The  noise  above  150°K  (close  to  the  Debye  temperature)  is  of  the  form  1/f  , 

and  requires  probably  a  vacancy  diffusion  model  or  other  specific  models  as  in 
the  work  of  Dutta  and  Horn.  Below  150°K  the  noise  is  1/f,  with  a  maximum  near 
80K,  and  then  a  continued  decrease.  If  the  N  in  Hooge's  formula  is  reinterpreted 
as  the  number  of  electrons  within  2  kT  from  the  Fermi  surface,  rather  than  the 
total  number  of  valence  electrons,  then  the  noise  can  well  be  understood  in  terms 
of  quantum  1/f  noise,  involving  Umklapp  processes. 

The  noise  in  GaAs  n+n  n+  mesas  of  submicron  dimensions  is  very  low.  The 
Hooge  parameter  is  of  order  10~^ ,  indicating  that  collisions  are  nearly  absent. 

In  l.lp,  n+n  n+  structures  intervalley  electron  transfer  is  noticeable.  An 
attempt  to  measure  the  associated  intervalley  noise  is  being  made.  The  n+p  n+ 
have  a  lot  of  noise  associated  with  the  prepunch-through  current.  This  is 
attributed  to  recombination  of  injected  electrons  via  empty  acceptors,  since  in 
the  unexcited  specimen  there  are  no  holes  due  to  electron  spillover.  At  very 
high  currents  (1-10  amps)  requiring  pulsed  noise  measurements,  ballistic  behavior 

with  low  noise  takes  over  again,  and  the  behavior  becomes  similar  to  that  in 

+  —  + 

n  n  n  structures. 


In  an  entirely  different  experiment  we  observed  for  the  first  time  a  1/f 


noise  component  in  radioactive  a-particle  decay  from  241  Americium.  This  noise 
was  deduced  from  counting  statistics  using  the  Allan  variance  theorem. 

Partition  1/f  noise  in  pentodes  was  measured  by  eliminating  the  cathode 
1/f  noise  with  feedback.  Van  der  Ziel  gave  an  expression  which  modifies  Handel's 
theory  for  this  type  of  noise.  The  observed  noise  seems  to  be  in  good  argument 
with  Handel's  modified  theory. 


ii 


In  the  theoretical  sections  some  further  remarks  were  made  on  the  inter¬ 


pretation  of  Hooge's  formula  for  metals.  Further,  we  calculated  quantitatively 
the  mobility-fluctuation  noise  associated  with  impurity  scattering  for  silicon 
and  for  gold.  It  is  found  to  be  extremely  low,  in  accordance  with  the  experi¬ 
mental  observations  that  impurity  scattering  gives  little  or  no  noise.  Computa¬ 
tions  for  acoustic  phonon  scattering,  involving  both  Normal  and  Umklapp  processes, 
are  being  carried  out.  Further,  a  survey  is  presented  of  the  particular  features 
of  nonequilibrium,  stationary  Markovian  noise  processes. 


iii 


1.  I-V  characteristic  for  n+p  n+  mesa;  p  layer  is  0.47p.  The  dashed  curves 
represent  the  results  computed  in  the  absence  of  recombination  current 
("prepunch- through"  current). 

2.  The  D.C.  conductance  for  the  same  sample  at  various  cryogenic  temperatures 
and  room  temperature. 

3.  The  noise  of  the  above  sample  at  300K. 

4.  The  noise  of  the  above  sample  at  77K. 

5.  Noise  of  a  0.24u  n  n  n  sample  at  300K. 

6.  Ibid,  at  77K. 

7.  Noise  of  a  0.4y  n+n  n+  sample  at  300K. 

8.  Noise  of  a  l.ly  n+n  n+  sample  at  300K. 

9.  Ibid,  at  77K. 
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10.  S^I  f  versus  current  (a  measure  for  Hooge's  a),  0.24y  sample. 

11.  Ibid,  for  0.4p  sample. 

12.  Ibid,  for  l.lVJ  sample. 

13.  Ibid,  vs.  widths  of  the  n  layers. 
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14.  I-V  characteristic  of  0.4u  n  n  n  sample. 

15.  I-V  characteristic  of  1.1  n+n  n+  sample;  note  the  intervalley  transfer. 

16.  Proposed  equipment  for  intervalley  high  frequency  noise  measurements. 

17.  Ibid. 

18.  Pulsed  noise  measurement  system  designed  by  Whiteside  and  Van  der  Ziel. 

19.  Pulsed  noise  results  for  n+p  n+  sample  (p~layer  of  0.47u). 
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20.  Sjl  f  at  4  MHz  for  n  p  n  sample  (pulsed  noise  measurement). 

21.  Ibid,  at  8  MHz. 

22.  Layered  structure  of  GaAs/GaAlAs. 

23.  I-V  characteristics  for  the  above  heterojunction  structure. 


24.  Noise  results  for  the  above  heterojunction  structure. 

25.  Counting  statistics  setup  for  a  particles  from  241  Am. 

26.  Energy  spectrum  for  a  particles  from  241  Am. 

27.  Convergence  of  relative  Allan  variance  with  increasing  number  N 
of  adjacent  counting  intervals;  T  =  1  minute. 

28.  Ibid;  T  -  3  minutes. 

29.  Relative  Allan  variance  versus  1/T  for  time  counting  intervals  up  to 
3000  minutes. 

30.  Allan  variance  versus  T  for  time  counting  intervals  up  to  3000  minutes. 

31.  Pentode  with  cathode  feedback  and  screen-grid  transformer  feedback. 

32.  Ibid,  in  more  expanded  fashion. 

33.  Current  noise  of  6  AU  6  WC  at  three  anode  currents. 

34.  Current  noise  of  6 AU  6 WC  ,  showing  effect  of  cathode  resistor  feedback 
(fixed  anode  current  of  3.62  mA) . 

35.  Ibid,  for  6  CEG  . 

36.  ST/ST  versus  R  for  6AU6WC  at  100  Hz. 
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37.  Ibid,  at  1  KHz. 

38.  Ibid,  for  6  CE  6  at  100  Hz. 

39.  Ibid,  at  1  KHz. 

40.  Same  results  for  6 AU  6 WC  in  log-log  scale  vs.  1  +  g  R  at  100  Hz. 
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41.  Ibid;  ordinate  is  now  i^/i  2 
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42.  Ibid  for  6  CE  6. 

43.  Current  noise  of  6 AU 6 WC  for  various  values  of  R  . 
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44.  Current  noise  of  6 AU 6 WC  for  two  different  anode  voltages. 

45.  Computed  relative  mobility-fluctuation  noise  due  to  impurity  scattering 
versus  T  in  silicon  films. 

46.  Similar  computation  for  gold  films. 
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INTRODUCTION 


In  our  previous  annual  report  (June  1982)  we  mentioned  six  types  of 
theories,  to  explain  1/f  noise,  that  have  survived.  Of  these,  we  can  now  omit 
the  category  "transport  noise  theories,"  a  specific  example  of  which  is  tem¬ 
perature  fluctuation  (or  heat-diffusion)  noise.  The  measurements  of  J.  Kilmer, 
performed  this  year  under  the  contract,  have  definitively  shown  that  the  noise 
in  electrically  isolated,  but  thermally  closely  connected>  metal  films  is 
uncorrelated  over  the  entire  range  300K-10K.  Other  transport  theories  had 
already  been  discarded  earlier  due  to  the  work  of  Mehta  and  Van  Vliet  (Physica 
Status  Solidi  (b)  106 ,  11,  1981).  So  far,  no  mathematical  transport  model  has 
yielded  a  1/f  spectrum  over  many  decades,  and  the  physical  evidence  for  the 
occurrence  of  transport  noise  as  a  basis  of  1/f  noise  is  entirely  lacking. 

Secondly,  in  the  June  1982  report  we  mentioned  the  model  based  on  phonon- 
distributed  lifetimes,  proposed  by  Jindal  and  van  der  Ziel.  So  far,  however, 
there  has  not  been  any  evidence  for  this  model  either,  while  the  occurrence  of 
very  long  phonon  lifetimes  remains  very  enigmatic.  So,  also  ruling  out  this 
theory,  four  types  of  theory  remain: 

1.  The  universal  theories; 

2.  The  van  der  Ziel-Bernamont-du  Pre-McWhorter  theories,  involving 
a  t  a(ci  w  1)  distribution  of  time  constants; 

3.  Specific  noise  model  theories; 

4.  The  mobility-fluctuation  bulk  model,  based  on  quantum  1/f  noise 
(Handel,  Ngai,  Widom,  et  al.). 

The  "universal  theories"  usually  seek  some  mathematical  general  model, 
which  explains  the  absence  of  time  scaling  (lack  of  characteristic  time  constants) 
in  1/f  noise.  For  example.  Montroll  and  Schlesinger  have  explained  that  a  log¬ 
normal  distribution  in  leading  order  distributes  a  random  variable  x  according 
to  x  ^  (Proc.  Natl.  Acad.  Sciences  USA  22.,  3380-3383,  1982).  Recently,  Marzec 


-  2  - 


and  Spiegel  offered  another  interesting  mathematical  example.  They  showed  that 
a  distribution  of  tunable  oscillators,  having  an  action  J  on  the  system,  and 


if  P^(w)  =  N/ui,  where  N  is  the  number  of  oscillators  (note  submitted  to  Nature, 
1983).  They  proceed  to  give  examples  of  the  occurrence  of  such  systems  of 
oscillators  in  physical  systems.  We  believe  at  present  that  such  universal 
theories  should  not  be  dismissed,  in  particular,  since  the  other  possibilities 
have  also  their  limitations. 

The  van  der  Ziel-  .....  -McWhorter  theories  still  are  viable  for  a  number 
of  cases.  In  these  theories  there  is  necessarily  always  a  lower  frequency 

=  1/t^  and  an  upper  limit  =  I/t^  outside  which  the  1/f  law  does  not  hold. 
Thus,  these  theories  are  not  scale-free.  In  some  cases  conclusive  evidence  for 
this  mechanism  has  been  found.  Thus,  Hanafi  and  van  der  Ziel  found  that  the  lower 
limit  could  be  varied  by  sputtering  off  layers  of  cadmium  mercury  telluride 
crystals.  We  note  that  the  noise  in  this  case  stems  from  the  surface,  as  in 
McWhorter's  original  experiments,  and  is  not  a  bulk  effect.  We  do  not  know  of 
any  physical  evidence  of  the  time-constant  distribution  model  in  cases  where  the 
noise  is  clearly  a  bulk  effect. 

The  "specific  noise  models"  certainly  apply  in  certain  cases.  Kilmer 

1  2 

found  that  the  noise  in  gold  films,  above  150K,  has  a  1/f  ’  spectrum;  the 
noise  rises  sharply  with  increasing  temperature  up  to  a  certain  maximum,  as  had 
also  been  observed  for  other  metals  by  the  Chicago  group  (Dutta,  Horn,  and 
others).  They  suggested  a  model  involving  vacancy  migration  for  this  temperature 
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range.  Likely,  some  effect  of  this  nature  is  responsible  for  the  noise.  Th *re 
are,  of  course,  a  host  of  other  "specific  noise  model"  theories,  such  as  the 
recombination  model  by  Min  (Solid  State  Electronics,  1979)  and  the  Island  model 
of  Pellegrini  (Phys.  Rev.,  1981).  However,  there  is  not,  in  our  opinion,  any 
evidence  for  the  validity  of  these  models.  And  if_  examples  for  the  occurrence  of 
such  theories  are  found,  then  still,  most  1/f  noise  cannot  be  explained  by  these 
models . 

The  quantum  theory  of  1/f  noise  as  advanced  by  Handel  in  1975,  and  worked 
out  in  more  detail  in  1980  (Phys.  Rev.  A22 ,  745,  1980),  remains  the  most  general 
theory  to  explain  1/f  noise.  However,  the  noise  calculated  from  Handel's  theory 
is  extremely  low.  Thus,  indeed,  this  theory  explains  at  most  the  "flicker  floor" 
occurring  in  physical  systems.  For  many  phenomena,  the  observed  noise  exceeds 
Handel's  theoretical  results  by  orders  of  magnitude.  In  this  report  we  present 
evidence  that  Handel's  theory  may  explain  1/f  noise  in  radioactive  decay,  in 
metal  films  below  the  Debye  temperature,  and  it  can  explain  partition  noise  in 
pentodes.  We  note,  however,  that  Handel's  work  still  receives  much  skepticism 
by  many  investigators.  It  has  to  be  admitted  that  a  number  of  points,  in  par¬ 
ticular  the  coherence  of  various  events,  responsible  for  the  1/N  dependence  in 
Hooge's  law,  is  not  well  understood  in  Handel's  theories.  This  makes  a  quantita¬ 
tive  comparison  with  observed  data  very  difficult,  as  we  will  see  in  this  report 
for  the  radioactive  decay  noise  and  for  the  1/f  partition  noise  in  pentodes. 

In  summary,  it  appears  to  us  that  the  explanation  of  1/f  noise  is  narrowed 
down  to  four  types  of  possible  theories.  The  time-constant  distribution  theory 
is  likely  only  applicable  to  surface-generated  noise.  The  specific  models  have, 
no  doubt,  validity  for  certain  types  of  noise,  in  particular  when  the  spectrum 
is  not  exactly  1/f,  but,  e.g.,  1/f®’^  -  1/f^"^.  Thus,  the  remaining  noise 
theories  of  promise  are  the  "universal  theories"  and  the  "quantum  theories"  of 
1/f  noise  which  involve  infrared  divergencies  due  to  nonlinear  coupling  of  the 
system  with  the  wave  field. 


A  point  which  also  needs  further  investigation  is  the  fact  that  order 
of  magnitude  variations  can  occur  for  macroscopically  identical  systems.  This 
was  abundantly  illustrated  at  the  Montpellier  International  Conference  (1983) 
for  the  noise  in  metal  films,  both  by  Fleetwood  and  Giordano,  and  by  our  own 
findings  due  to  Kilmer  et  al.  Why  seemingly  "identical"  systems  can  have  dif¬ 
ferent  magnitude  of  noise  (a  factor  3-10)  is  at  present  an  enigma. 

B.  EXPERIMENTAL  WORK 

Much  experimental  work  was  presented  at  the  Seventh  International 
Symposium  on  the  Noise  in  Physical  Systems  and  Third  International  Conference 
on  1/f  Noise,  held  in  Montpellier,  France,  May  17  -  20.  Where  appropriate,  we 
reproduce  the  contributions  of  this  conference. 


I.  1/f  NOISE  IN  METAL  FILMS  OF  SUBMICRON  DIMENSIONS 


J#  Kilmer,  C.M.  Van  Vliet,  G.  Bosman  and  A.  van  der  Ziel 


Department  of  Electrical  Engineering 
University  of  Florida,  Gainesville,  FL  32611,  USA 


Photonic  quantum  1/f  noise  has  been  identified  in  Au  metal  films  below  the  Debye 
temperature.  The  low  values  of  the  Hooge  parameter  predicted  by  Handel's  theory 
(i.e.,  *true  mr  10”6  to  10“ 8 )  are  arrived  at  by  realizing  that  only  the  fraction 
3kT/2L£  y  of  total  number  of  carriers  are  available  for  scattering  at  a  given 
temperature. 


In  recent  years  considerable  progress  has  been 
made  in  the  understanding  of  1/f  noise.  It  is 
now  well  established  that  in  many  cases  there 
is  fundamental  1/f  noise  caused  by  mobility 
fluctuations,  in  particular  by  fluctuations  in 
the  scattering  cross  section  of  scattering  of 
electrons  by  phonons. [1) [2} 


The  only  general  theory  of  1/f  noise  which  can 
explain  such  fluctuations  was  given  by  Handel 
in  1975.(3] [4],  However,  until  recently  experi¬ 
mental  evidence  verifying  the  theory  did  not 
exist.  Specifically,  Handel’s  quantum  i/f 
noise  theory  was  questioned  as  the  source  oi 
1/f  noise  in  electronic  circuits  because  of  the 
low  value  of  the  Hooge  parameter,  ajj,  calcula¬ 
ted  from  his  theory.  Brieily,  the  theorv  states 
that  the  interference  between  the  part  of  the 
carrier's  wave  function  which  suffers  lasses 
due  to  an  inelastic  or  "brensst rahlung"  scat¬ 
tering  under  the  emission  of  infraquanta  and 
the  part  of  the  wave  function  which  does  not 
suffer  losses  produces  very  low  energy  beats 
which  translate  themselves  (At  = hf )  as  1/f 
noise.  Handel's  theory  predicts  that  scatter¬ 
ing  involving  Umklapp  processes  (U-processes) 
provides  the  largest  source  of  1/f  noise  in 
metals  since  the  photon  infraquanta  coupling 
constant,  aA,  is  given  by[M 


lA 


(1) 


where  a  is  the  fine  structure  constant  (137)” ‘, 
perhaps  modified  by  the  dielectric  constant  in 
the  metal,  c^is  the  velocity  of  light  in  the 
metal,  and  \p/fi  3  '.k  is  the  change  in  wave  vec¬ 
tor.  Since  the  U-process  gives  the  largest  Ik, 
we  expect  them  to  be  the  largest  contributor  to 
1/f  noise  in  metals.  Though  the  dielectric 
constant  of. metals  is  not  well  known,  and  may 
he  complex,  one  e.is i  ly  sees  that  the  correc¬ 
tions  in  »  and  e  cancel,  so  that  we  can  turthcr 
take  the  free  space  values. 


Gold  thin- film  resistors  (2,000  A  thick)  were 
prepared  for  us  by  Dr.  E.  Noll  and  R.A.  Buhrman 
of  the  National  Research  and  Resource  Facility 
for  Submicron  Structures  at  Cornell  University. 


fhe  length  of  the  resistors  is  close  to  800  urn, 
and  we  have  measured  the  1/f  noise  in  1  ym-width 
samples.  These  dimensions  give  a  resistance  of 
a  few  hundred  ohms  and  the  noise  spectrum  can  be 
readily  measured,  after  amplification,  by  an 
HP  3582  Spectrum  Analyzer.  By  incorporating  a 
calibrated  noise  source,  the  absolute  magnitude 
of  the  resistor's  current  noise  spectrum,  St, 
can  be  directly  calculated  by  simply  comparing 
the  relative  spectra  of  device  on,  device  off, 
and  calibration  source  on;  it  can  be  shown  that 
the  amplifier's  parameters  cancel  out.  The  gold 
films  were  mounted  to  the  cold  head  of  a  CTI 
Cryogenics  Model  21  liquid  He  closed-cycle  re¬ 
frigerator  capable  of  maintaining  a  stable  tem¬ 
perature  (i.e.,  -‘0.1  K  over  the  duration  of  a 
low- j rcquency  noise  measurement)  anywhere  be¬ 
tween  300  K  and  3  K. 

The  results  of  the  experiment  give  current  spec¬ 
tra  proportional  to  1/f  between  1  and  100  Hz. 
Below  L  Hz  we  have  cryostat  noise,  and  above 
100  Hz  wo  have  device  noise  competing  with  the 
amplifier's  noise  and  thermal  noise.  The  slope 
shows  y  sr  1.2  from  300  K  to  about  the  Debvc  tem¬ 
perature  (0q  =  165 K  for  gold)  in  agreement  with 
Fleetwood  and  Giordano. (5]  Belov  the  Debye  tem¬ 
perature  the  slope  evens  off  to  y  —  1.  indicat¬ 
ing  a  more  "pure"  1/f  noise  present  at  the  lower 
temperatures.  Next,  we  characterize  the  magni¬ 
tude  of  the  1/f  noise  by  calculating  the  dimen¬ 
sionless  Hooge  parameter,  according  to  the 

formula 

s  (f)  a 

_ L _ -  _  f  ”>\ 


where  N  is  the  total  number  ot  available  elec¬ 
trons  in  the  metal;  for  those  spectra  where  . 
was  *  1,  we  took  f  =  10  hz.  In  Figure  1  ve  have 
used  N  =  nV  where  we  set  n  ss  IU”'  cm”  and  V  is 
the  sample  volume  (300  ;.m  *  1  umx0.2  wn). 

Figure  1  shows  an  interesting  dependence  of  the 
magnitude  ol  below  the  Debye  temperature 
corrosponu ing  to  the  occurrence  of  the  more 
"pure”  l/f  noise. 


10 


,-2J 


kT3(log  N)/3£p  where  £  p  is  the  Fermi  energy. 
Explicitly ,  we  have 


kT9(log  N)  _  Jji. 1  (  ~  ^F)/kT  1 

3eF  ’  T.i,  [(<£<.- <5F)/k T] 


! 

I 

I 


/O  5  L.. 
'OK 


!OCK  f  300K 
Qo 


where  £  is  the  bottom  of  the  conduction  band, 
and  k  is  the  Fermi  integral  of  order  k.  For 
total  degeneracy,  ^^(n)  =  n*  +  Vr(k  +  2)  .  Thus 
the  ratio  (4)  becomes  2A£F/3kT  with 
A  ^p  =  £p-  £c.  Consequently,  eq.  (3)  followed 
by  the  proper  statistical  summation  leads  to 


sx(f) 


i!£_F 

3kT 


fN 


(5) 


indicating  that  the ^number  of  carriers  available 
for  scattering  is  N"  *  N( 3kT/2A£  p) .  This  is 
also  intuitively  obvious:  the  Fermi  function 
differs  only  appreciably  from  1  or  0  in  a  slice 
of  order  kT.  That  such  a  reduction  in  noise 
must  occur  in  metals  was  perhaps  first  pointed 
out  in  a  classic  paper  by  Brillouin[ 111  on  the 
first  noise  observations  in  metals,  by 
Bernamont . [ 12 1  Comparing  now  (5)  with  (3)  we 
find  that  the  "true"  Hooge  parameter  is  related 
to  the  observed  Hooge  parameter  by 


=  3kT 

atrue  2A£  p  *H* 


(6) 


Figure  1.  a  as  a  function  of  temperature 


With  1  £,  p  =  5.5  eV,(13]  the  values  of  -«true 
were  computed  to  yield  the  data  of  Fig.  J . 


We  believe,  however,  that  (2)  is  not  a  proper 
characterization  of  noise  in  metals.  As  was 
poincc  1  out  by  Van  Viiet  and  Zijlstra,[6)  the 
basic  formula  for  the  mobility  fluctuations  for 
th  scattering  of  a  single  carrier  is 


yf) 

u/ 


true 

f 


(3) 


where  we  subscripted  the  j-value  as  a  .To 
r  true 

obtain  the  fluctuations  in  the  band  nobility, 
or  for  that  matter  of  the  current  I,  we  must 
sum  over  the  scattering  fluctuations  of  all  car¬ 
riers  in  the  band. [7] [3]  For  a  nondegenerate 
semiconductor  this  provides  a  factor  1/N  in  the 
denominator,  see  12).  In  metals,  however,  most 
of  the  carriers  are  "frozen”  in  the  Form  sea, 
a  fact  also  noted  by  Dutta  and  Horn  in  their 
review  paper. [9]  Therefore,  as  we  pointed  out 
elsewhere,  in  connection  with  thermal  noise  from 
metals  and  the  Einstein  relation, [ 10 1  the  result 
of  the  summation  must  lie  multiplied  by 
<  AN’  >  /  <N>  «  where  the  subscript  "  "  refers 
to  the  grand  canonical  ensemble. | 10 ,  sec.  1] 

From  statistical  mechanics  the  above  factor  is 


As  is  noted,  we  now  obtain  -^-values  low  enougn 
to  become  in  the  ballpark  expected  from  the 
quantum  theory  of  1 / f  noise.  In  the  latter 
theorv,  the  '*Crue  of  eq.  (3)  is  just  twice  the 
infrared  exponent,  i.e.. 


atrue 


2a  A. 


(7) 


With  Ap/m  -  2vpsin  t/2,  where  Vp  is  the  Fermi 
velocity  (1.39*  10-3  cm/sec)  and  ?  is  the  scat¬ 
tering  angle  (~150°  for  L’-procosses) ,  we  ob¬ 
tain  from  (1)  2aA  s:  2.4 x  10“  .  This  value  is 
approximate  since  N  is  not  exactly  known  and 
since  more  correctly  we  must  take  into  account 
the  detailed  geometry  of  the  Fermi  surface,  be¬ 
ing  a  sphere  with  eight  "necKs"  (see  ZinanJlA]) 
However,  this  value  comes  close  to  the  observed 
valuo  of  (.»tril0)max  in  Fit;.  2,  being  '..9^10'  . 

Qual itativelv ,  we  believe  that  the  observed 
data  of  Fig.  2  can  be  well  understoou.  Above 
the  Debve  temperature  y,  region  C,  sene  non- 
fundamental  l/f'  noise  occurs,  similar  to  the 
"tvpe  b"  noise  observed  bv  Dutta  and  Horn. [9) 
Below  we  have  ror  the  first  time  a  clear 
indication  of  the  occurrence  or  quantum  l/f 
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Figure  - .  i  as  a  function  of  temperature 

true 


noise.  (.’-processes  dominate  in  the  region  B, 

In  the  region  A,  U-processes  freeze  out  and 
normal  phonon  processes  (N-processes)  take  over. 
Finailv,  in  a  region  i)  (not  yet  observed), 
ionized  impurity  scattering  nay  give  rise  to  a 
floor  at  very  low  temperatures# 


A  quantitative  theory  has  not  yet  been  fully 
developed.  However,  with  (.’-processes  dominat¬ 
ing  the  noise, and  N-processes  and  ionized  in- 
pur i tv  scattering  dominating  the  resistance, 
one  expects  the  temperature  dependence  to  bo  of 
t ae  torn 


wnere  ..  i:.  the  i  nki.ipp  tt  npor  lturo, 

*  hq  a  •  ,  ,  in  whicn  a  is  the  transverse 
velocity  'i  =ouna  in  gild,  ij  is  :*••!  i  zn.imi  ’  •: 

onst.uit,  r.d  •;  ;  .  .»  pisonon  vector  associate*! 

’.:ith  t  tie  iengtn  i  the  no*  ka  tween  id  a. ant 
.Vrni  an  i.  >  i:\  tin-  extended  c.  ne  scl.«r.e,  .*e 

onputeJ  a  ■;  *i.j  I'l  .  le 


uQ  =  1.2  *  10 5  cm/ sec.  This  yields  Gji  *  S7.6K. 
The  observed  maximum  in  the  noise  occurs  at 
about  »>0  K .  Though  nanv  details  need  fuller 
consideration,  we  believe  that  the  observed 
noise  can  be  reasonably  well  explained  by  the 
proposed  processes.  Development  of  a  full 
theorv  and  detailed  1,'f  noise  measurements  in 
metals  nuv  give  much  insight  into  the  nature  of 
the  phonon  processes  undergone  by  the  Fermi 
surface  electrons. 

We  finally  note  that  another  type  of  confirma¬ 
tion  of  the  theory  of  quantum  1/f  noise  was 
recently  provided  by  a-particio  decay  statis¬ 
tics.  [  15  ] 
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Measurements  on  the  0.5  pm  and  2  pm  width  gold  thin  film  samples 


are  presently  being  performed  and  a  similar  a  vs.  T  plot  is  observed. 

Also,  to  get  a  more  complete  picture,  the  same  noise  measurements  as  a 
function  of  temperature  will  be  performed  on  other  metal  films  with  different 
Debye  temperatures.  Presently,  masks  are  being  prepared  to  evaporate  aluminum 
on  a  silicon  oxide  substrate  in  a  5  urn  *  1000  pm  *  2000X  thin  film  resistor 
configuration  here  at  University  of  Florida's  microelectronics  laboratory. 
Producing  films  of  less  than  5  pm  width  will  be  aimed  for,  since  we  desire  a 
greater  resistance  from  the  thin  film  resistor  samples.  Aluminum  is  not  ideal 
since  its  0^  is  428 K  and  the  resistor  must  be  placed  in  an  oven  to  probe  the 
noise  near  the  Debye  temperature. 

After  perfection  of  the  aluminum  evaporation,  palladium  (9  =  274  K)  , 

silver  (9  =  225  K)  ,  chromium,  and  indium  which  can  be  evaporated  will  be  used 

to  produce  resistors  in  the  same  thin  film  configuration,  and  noise  measure¬ 
ments  will  be  performed. 

Hopefully,  we  will  be  able  to  produce  platinum  films  (0^  =  240  K)  in 
our  microelectronics  lab  when  we  get  the  new  ion-beam  sputtering  equipment. 

Finally,  provisions  have  been  made  to  make  noise  measurements  on  samples 
in  the  8K  to  2 K  ambient  temperature  range  by  directly  submerging  the  sample 
in  a  liquid  helium  cryostat  and  reducing  the  vapor  pressure  by  pumping. 
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II.  NOISE  IN  NEAR-BALLISTIC  n'nn"  AND  n' pn* 
GALLIUM  ARSENIDE  SUBMICRON  DIODES 

R.  R.  Schmidt.  G.  Bosman  and  C.  M.  Van  Vlilt 
Department  of  Electrical  Engineering.  University  of  Florida.  Gainesville.  FL  32611  U  S  A 

and 

L.  F.  Eastman  and  M  Hollis 

School  of  Electrical  Engineering.  Cornell  University.  Ithaca.  NY  14853.  L'.S.A 

{Receited  17  August  I ‘>82.  in  reused  form  26  October  1982) 

Abstract  Direct-current  Id  c  I  characteristics  and  noise  measurements  in  the  range  I  Hz  2 5  k H /  are 
reported  for  n  '  mi  '  and  n  '  pn  '  near-ballistic  devices,  with  n  regions  tp  regions)  of  0.4  jim  (0.45  /im|. 
fabricated  by  molecular  beam  epitaxy  at  Cornell  The  n  '  nn  ‘  mesa  structures  show  very  low  I  /  noise, 
indicating  a  Hooge  parameter  ia  =  6!)«  10  ’  This  very  low  noise  is  attnbuted  to  the  near  absence  of 
phonon  collisions  The  thermal  (  -  like)  noise  above  1  kHz  is  equal  to  Nyquisl  noise  at  the  lowest  currents, 
rising  to  slightly  above  Nyquist  noise  for  high  currents,  indicating  the  presence  of  carrier  drag  effects.  The 


n  ’  pn  '  noise,  on  the  contrary,  is  quite  high  It  seems 
for  low  injection  of  electrons  in  the  p  region.  The 
ballistic  or  near-ballistt :  behavior  is  discussed 

I.  INTROfH  CTION 

Submicron  gallium  arsenide  structures  are  of  great 
current  interest,  since  they  permit  ballistic  or  near- 
ballistic  electron  flow,  which  in  turn  leads  to  carrier 
velocities  that  far  exceed  the  saturation  velocity  in 
collision-dominated  conduction,  thus  enabling  the  de¬ 
sign  of  picosecond  sw  itching  devices  and  other  novel 
applications.  The  fabrication  of  submicron  devices 
has  been  made  possible  by  modern  MBF.  techniques, 
electron  lithography,  etc.  For  GaAs  near-ballistic  be¬ 
havior  requires  that  the  distance  to  be  traveled  by  the 
injected  electrons  is  less  than  or  of  the  order  of  0.7  ft  m. 
Eastman  ct  til.  report [  I  j  that  the  mean  free  path  for 
phonon  emission  into  optical  polar  modes  at  room 
temperature  is  0.1  jim  for  electrons  of  0.05  eV.  and 
0,2  qm  for  electrons  of  0  5  cV.  Phonon  absorption  has 
a  longer  mean  free  path  and  can  be  neglected  for  the 
devices  reported  here,  having  thicknesses  of  0  4/rm 
(thickness  of  n  layer  in  n  '  nn  '  devices)  and  0  45  p m 
(thickness  of  p  layer  in  n  '  pn  '  devices).  At  higher 
electron  energy  intervalley  scattering  becomes  im¬ 
portant.  thus  limiting  the  near-ballistic  range  to  about 
0.5  cV  of  electron  energy  In  a  sample  of  0.4  jim  thick¬ 
ness  about  two  phonon  emissions  may  occur.  These 
involve,  however,  small  angle  deflections  only  (5-10  ) 
and  have  little  effect  on  the  d.c.  carrier  characteristics, 
according  to  Ref  [I] 

The  theory  for  "pure"  ballistic  behavior  (no  col¬ 
lisions  sulfered  whatsoever)  was  developed  by  Shur 
and  Eastman  in  1979  m  a  basic  paper  on  this  topic [2). 
Thev  solve  Poisson's  equation,  allowing  for  space 
charge ot  both  lived  tont/ed  donors(or  acceptors)  and 
iniccted  carriers  Employing  boundary  conditions 
which  neglect  the  initial  thermal  energy  of  the  elec¬ 
trons  iniccted  Irom  the  n  '  into  the  n  layer  (or  p  layer 
alter  punch-through l.  they  lind  the  solid  state  analog 


to  be  associated  with  the  ambipolar  effects  occurring 
importance  of  noise  measurements  for  confirming 

of  Child's  law  in  vacuum  tubes.  For  sufficiently  small 
voltages  there  is  a  domain  in  which  the  current  /  goes 
as  I'1 when  the  injected  space  charge  exceeds  the 
fixed  charge  due  to  the  ionized  donors  or  acceptors, 
the  characteristic  changes,  however,  to  the  familiar 
r'-  form  In  Ref. 1 1 )  measurements  are  presented 
which  fairly  well  support  these  predictions,  providing 
the  nonparabolicity  of  the  bands  and  the  onset  of 
intervalley  scattering  at  higher  voltages  are  taken  into 
account.  In  a  later  theory.  Shur(3]  and  Shur  and 
Eastman  [4]  extended  the  theory  to  that  for  "near- 
ballistic"  devices,  in  which  few  collisions  can  occur. 
Since  the  Boltzmann  equation  would  be  inappropriate 
for  that  regime,  the  collisions  are  taken  into  account 
by  adding  momentum  and  energy  "drag  terms "  to  the 
otherwise  ballistic  equations  of  motion.  In  this  way. 
the  transition  from  Child's  law  (r— >  x .  where  r  is  the 
collision  time)  to  the  Mott  and  Gurney  law.  /xl  ‘ 
(finite  r)  is  covered  by  this  approach. 

Two  modifications  have  been  proposed  by  others, 
which  may  have  a  bearing  on  the  present  paper.  First. 
Rosenberg  ct  u/.[5]  discuss  the  effects  of  "spillov  er"  of 
carriers  at  the  n  "  n  high-low  junction.  This  means,  in 
essence,  that  the  boundary  conditions  must  be 
changed  to  account  for  the  depletion  of  n  *  regions 
and  spillover  into  the  adjacent  Debye  lengths.  As  a 
result,  the  effective  width  of  the  n  region  is  smaller  and 
the  current  is  higher  than  that  computed  in  Ref.  [2] 
Secondly  .  Cook  and  Jeffrey  |6)  have  indicated  that  the 
energy  or  velocity  distribution  of  the  electrons  cannot 
be  neglected.  The  velocity  dispersion  is  accounted  for 
by  the  introduction  of  an  electron  temperature  gra¬ 
dient  term  in  the  momentum  balance  equation  {op  at 
eqn  (8)).  Though  they  argue  that  this  leads  to  the 
occurrence  of  a  potential  minimum  somewhere  be¬ 
yond  the  "cathode"— rather  than  at  the  cathode— 
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similar  to  Langmuir's  treatment  of  vacuum  diodes,  we 
have  great  reservations  about  their  treatment.  A  more 
straightforward  and  correct  approach  was  very  re¬ 
cently  presented  hv  Holden  and  Dehnev  [7],  in  a  paper 
based  on  ideas  from  the  well-known  vacuum  diode 
discussion  of  Fry  [N|  Their  results  indicate  that  no 
lived  power  law  can  be  stated  For  0.5 /j  samples  the 
limiting  slope  is  found  to  be  1.14.  For  the  noisecharae- 
teristics  the  inclusion  of  the  velocity  dispersion  is  es¬ 
sential.  in  particular  for  the  high-frequene>  thermal- 
like  velocitv -tluctualion  noise 

In  this  paper  we  describe  verv  accurate  low- 
frequency  and  high-frequency  noise  measurements  on 
near-ballistic  devices  Such  measurements  serve  a 
threefold  purpose  First,  from  a  practical  point  of  view 
noise  data  reveal  thepractic.il  performance  limitations 
of  the  novel  high-speed  devices  As  we  will  indicate 
the  noise  of  the  n  '  nn  '  devices  is  extremely  low.  the 
n  '  pn  '  devices,  however,  fare  much  worse.  Secondly, 
noise  measurements  at  audio  and  suhaudio  fre¬ 
quencies  shed  much  light  on  the  I  /  nt  ise  problem. 
According  to  most  recent  theories,  such  noise  is 
thought  to  be  caused  by  mobility  fluctuations  (see.  e.g. 
Hoogccr  al  [U]  and  van  der  /iel(IO|).  If  collisions  in  the 
near-ballistic  regime  are  rare,  one  expects  the  I  f  noise 
to  he  very  low  and  ultimately,  in  "pure"  ballistic 
devices,  to  he  absent  Our  work  on  n  '  nn  '  devices 
indicates  that  this  could  be  correct  Third,  and  not 
least,  we  believe  that  the  high-frequency  noise  (ther¬ 
mal.  velocitv -fluctuation,  or  diffusion  noise)  will  shed 
much  light  on  the  mode  of  operation  of  near-ballistic 
devices.  To  date,  no  full-fledged  theory  for  sue'  use 
exists;  we  only  have  some  preliminary  computations 
by  van  der  Ziel  and  Bosnian] 1 1.  12]  However,  once 
this  noise  is  understood,  we  will  have  a  powerful 
means  of  substantiating  or  amending  the  various 
theories  on  near-ballistic  behavior. 


m 


fig  I  /’  type  near-ballistic  mesa  structure; 

structure  is  similar. 


the  n-type 


2.  f  \pmi\iF\T  \i 

The  near-ballistic  diode  (NBD)  is  a  sandwiched 
mesa  structure  of  five  lightly  doped  p  or  n  layers, 
alternating  with  heavily  doped  n  '  layers,  sec  Fig  I 
The  doping  densities  of  the  various  regions  are 
10  'em  for  the  n  '  regains,  approx  2  *  10'"  cm 
for  the  n  regions  and  approx  10''  cm  '  for  the  p 
regions  The  diameter  of  the  mesas  is  100  qm  The 
devices  were  manufactured  by  molecular  beam  epi¬ 
taxy  at  the  Cornell  I  niversity  Submicron  Research 
Facility  The  mesas  were  provided  with  very  low 
ohmic  Xu  Ge  contacts.  A  low  frequency  equivalent 
circuit  of  the  n  '  nn  ‘  device  is  given  in  Fig  2  The 
main  element  <ri  regions)  has  a  resistance  of  0  75  Q 
To  infer  the  correct  noise  of  the  device,  the  values  of 
the  parasitic  resistances  in  Fig  2  have  been  taken  into 
account.  For  the  n  '  pn  '  devices  the  p  layer  gave  a 
resistance  of  order  00  Q  at  I  m.A.  the  parasitic  re¬ 
sistances  in  this  case  were  negligible 

The  characteristics  of  the  two  types  of  devices  arc 
quite  different  The  noise  measurement  of  the  n 
NBD's.  in  particular,  was  a  challenge  To  do  this  we 
used  the  setup  shown  in  Fig.  3  A  Hewlett  Packard 


lie  2  Fqimalcm  circuit  of  n  nn  '  structure,  showing 
parasitic  elements  The  contacts  (2)  and  @  refer  to  the  top 
contacts  I  in  fie  I.  contact  (j  refers  to  I  .  and  contact 
©  to  CiNO  of  l  1C  *  l 


He  J  Correlation  measurement  setup  **Ol  T“  is  device 
under  test  and  L\*\  denotes  the  two  linear  amplifiers  ol  the 
measurement  channels 


PAMl". 


11 


Noise  in  neur-ballisiic  gallium  diodes  J  W 


t 


Fig  4  I  V  characteristic  n  '  «n  '  device. 


F5S2  spectrum  analyzer.  featuring  a  dual-channel  last- 
Fourier  transform  method,  was  employed  B\  mea¬ 
suring  the  coherence  (square  of  the  correlation)  be¬ 
tween  the  two  channels,  noise  levels  significant!) 
below  the  noise  level  of  the  preamplifier  could  he 
detected  For  the  preamplifiers  we  used  live  common 
emitter  transistors  GE;S2  in  parallel  This  resulted  in  a 
7  il  noise  resistance  lor  frequencies  above  70  H/.  The 
equivalent  noise  resistance  of  the  cross-correlation 
setup  was  found  to  be  as  low  as  0  2  ft.  thus  enabling  us 
to  accurately  measure  the  thermal  noise  of  the  very 
low  ohnne  n  '  nn  ’  devices 
The  d.c.  /  I  characteristic  of  an  n-type  device  is 
shown  in  Fig  4  For  the  higher  voltage  low  duty-cycle 
pulsed  measurements  were  made  The  highest  voltage 
over  each  layer  is  about  0.2  V.  well  below  the  occur¬ 
rence  of  inters  alley  scattering  We  notice  that  w  ithin 
the  experimental  errors,  the  characteristic  is  entirely 
linear,  in  agreement  with  the  theoretical  curves  given 
in  Ref  |7]  (see  their  Fig.  2.  L  —  0.5 n  sample  with 
F  <  O  F  v  In  accord  with  their  conclusions,  we  note 
that  near-ballistic  behavior  cannot  usually  be  deduced 
from  just  a  measurement  of  the  /  F  characteristic. 
The  resistance  value  deduced  from  Fig  4  is  0.75  fl. 
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Fig  ?  I  V  characteristic  of  n  '  pn  '  device 
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Fig  6  Noise  spectra  for  n  '  nn  '  device. 

Alternating  current  (u.c.)  impedance  measurements  j.  noise  of  •*«•  -  device 

between  1  Hz  and  I  kHz  showed  a  flat  response  and  The  magnitude  of  the  noise  current  spectrum  for 
confirmed  this  resistance  value.  These  results  are  four  different  currents,  in  the  frequency  range 
similar  to  the  latest  results  reported  by  Hollis  i  t  at.  1 1 3).  1  Hz  25  kHz.  is  show  n  in  Fig.  6  Thermal  levels  and 

The  d.e.  /  C  characteristic  of  a  p-type  device  is  excess  1  /  noise  are  seen  To  determine  the  thermal 

shown  in  Fig.  5.  The  device  is  linear  up  to  a  current  (-like)  noise  levels,  the  1  /'components  are  subtracted, 

level  of  1  mA.  corresponding  to  about  /?,  =  90  fl.  The  The  results  are  shown  in  Fig.  7.  The  levels  are 

slope  then  increases  to  a  value  of  about  3  in  the  10-  averaged  over  the  frequencies  for  which  there  is  a 

100  mA  range.  Finally,  at  very  high  currents  the  slope  plateau  (1-25  kHz;  the  75  mA  curve  may.  however, 
becomes  less,  perhaps  approaching  three-halves  as  the  show  some  g-r  noise  from  1  7  kHz:  the  thermal-like 

slope  falls  olf.  The  a.c.  resistances  are  again  fiat  for  all  noise  occurs  for  7  kHz  and  higher).  The  ratio  of  these 

measured  frequencies  (up  to  100  kHz)  averages  io  4  kT  0.75  is  plotted  versus  bias  current  in 
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I  ig  ?  Thermal  Hike*  noise  for  n  *  nn  *  device 
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Fig  H  Plot  of  thermal  < -like >  noise  of  n  ‘  nn  '  diode  vs  diode  currcm. 


Fig.  8.  We  note  thut  there  is  an  indication  that  the 
noise  exceeds  the  true  thermal  noise  4 kT  R,  at  the 
higher  bias  currents. 

The  1  /  slope  of  the  noise  for  the  higher  current 
levels  is  clearly  seen,  and  straight-line  approximations 
arc  made  to  the  data.  The  values  so  obtained  at  10  hz 
are  plotted  versus  bias  current  in  Fig  9.  We  note  that 
the  expected  behavior  for  1  /  noise.  S,  s  /'.  is  well 
satisfied. 

4.  NOISF.  OF  n'pi i*  DEVICE 

The  noise  current  spectrum  versus  frequency  for 
several  bias  currents  of  an  n  '  p  n  ’  device  is  shown 
in  Fig.  10.  The  excess  low-frequency  noise  of  this 
device  is  orders  of  magnitude  larger  than  for  the 
n-type  device.  Another  notable  feature  is  the  fre¬ 
quency  dependence,  which  shows  a  slope  of/  °’  to 
f  "\  Extrapolating  to  the  corner  frequency  above 
w  hich  thermal  noise  dominates  gives  a  value  of  over 
100  MHz  for  even  the  lowest  (I00p.4)  bias  current. 

The  dependence  of  the  noise  current  on  bias  cur¬ 
rent  at  100  Hz  is  displayed  in  Fig.  1 1.  There  is  an  l: 
dependence  up  to  about  I  mA.  At  higher  currents  the 
noise  increases  less  fast  and  probably  goes  through  a 
maximum. 

S.  DISCISSION  OF  »•»*•  RF.SCI.TS 
(a)  I  /  now 

In  I'i(i9  llooge  developed  the  following  empirical 
formula  for  I  /  noise. 

SV(D  !'  =  t„  /.V.  (5.1) 

where  /  is  the  frequency.  ,V  the  total  number  of 


carriers  in  the  sample  contributing  to  the  noise,  and 
is  Hooge  s  parameter.  Initially.  iH  was  thought  to 
be  a  constant,  of  order  2  x  10  '.  Later  on.  it  was 
found  that  material  variations  for  i„  do  occur, 
w  hereas  in  addition  decreases  as  in  nX  if  impurity 
scattering  dominates  over  lattice  scattering  (/a,.): 
Bosman  cl  al  [14]  also  found  that  aH  decreases  due 


10  100 
I(m&) 

Fig  9  I  /  noise  of  n  ‘  nn  •  diode  vs  current  at  /  *  10  H: 
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f  ig  10  Noise  spectra  for  n  '  pn  '  device. 
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to  carrier  heating  In  a  nonhomogeneous  sample  in 
which  the  carrier  density  is  a  function  of  position. 
n(\).  such  as  occurs  in  our  mesas  due  to  spillover 
(Section  ll  and  injection,  eqn  1 5.1)  must  be  modified 
It  is  easily  shown  (van  der  Ziel  and  van  Vliet [ 1 5] I  that 
in  that  case  (5.1)  is  to  be  replaced  by 


4\/t t)  I'  -  (Uw  1.4 /. 


.  r  dv 

.  m)  - -. 

!„ «( v ) 


(5.2) 


where  m  is  the  number  of  n  layers  in  scries.  L  the 
width  of  one  >1  layer,  and  .4  the  cross  section.  This 
formula  is  correct,  whether  or  not  the  motion  of  the 
carriers  is  ballistic  Whereas  the  detailed  profile  n(  v  1 
is  complex  due  to  spillover  and  carrier  flow,  we  may 
assume  that  for  most  of  the  layer  rt(  v  t  >  1.5n„.  where 
n„  is  the  doping  dcmsty|l6)  With  this  estimate  we 
obtain  for  iH  from  Fig  Id.  i„  =  6  *  10  '  This  is 
considerably  lower  than  the  value  of  i„  for  bulk 
CiaAs,  which  Hoogc  or  u/  |P]  list  as  6  »  10  ’  Thus 
the  measurements  of  this  paper  confirm  that  col¬ 
lisions  arc  mainly  absent  in  this  device  Moreover,  if 
Handel's  theory  of  I  f  noise  is  vahd|ls|.  verv  low 
noise  can  be  expected  from  those  collisions  which  still 
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occur.  involving  polar  phonon  emission.  As  we 
noticed.  the  defection  angle  (I  for  such  processes  is 
\er>  small,  whereas  in  Handel's  theory  of  quantum 
1  /  noise  the  magnitude  goes  as  sin;  :0.  Measurements 
on  0.14pm  devices  arc  underway  Very  recent  data 
by  J  Andrian  < 7  <i!  [  1  d)  indicate  that  for  these  devices 
t,i  shows  a  continued  decrease.  So  far.  these  results 
are  the  best  confirmation  yet  that  I  /  noise  is  caused 
by  lattice  phonon  collisions. 

(b)  Thermal  noise 

The  designation  "thermal  noise"  is  used  here  for 
the  thermal-like  noise  observed  at  high  frequencies. 
In  a  collision-limited  device  this  noise  is  due  to  the 
dilTusion-noise  source,  which  by  Einstein's  relation 
transforms  to  a  thermal-noise  source  for  cold  elec¬ 
trons.  In  the  space-charge  limited  injection  operation 
(Mott  Gurney  law),  the  noise  becomes  then  8kT  /?,. 
[see  Ref.  Ill  In  a  pure  ballistic  device,  on  the  other 
hand,  this  noise  is  due  to  shot  noise.  However,  the 
vacuum  case  shows  that  the  noise  is  distinctly 
governed  by  the  velocity  distribution  of  the  emitted 
particles  Thus,  a  treatment  as  the  Child's  law  analog 
of  Ref.  1  will  not  suffice  to  obtain  the  noise:  the  latter 
must  be  patterned  after  Langmuir's  derivation  of  the 
d.c.  characteristic;  sec  in  particular  the  noise  treat¬ 
ment  by  D.  O.  North [10]  and  Schottky  and 
Spenke[!l).  Lacking  a  detailed  theory,  van  der  Ziel 
and  Bosmanjll]  indicated,  nevertheless,  that  sub- 
thermal  noise.  (W  T  R,  with  0  <  \.  can  be  expected. 
This  is  not  corroborated  by  the  results  of  Fig.  8. 
While  it  is  very  unlikely  that  the  collision-limited  case 
applies  in  view  of  the  low  |  /  noise  reported 
above  it  is  likely  that  carrier  drag  effects,  such  as 
considered  in  Refs.  (.1.  4.  b],  take  place.  These 
effects  should  be  incorporated  by  considering  a  Lan- 
gevin  equation  patterned  after  the  momentum  and 
energy  balance  equations  of  Shur  and  Eastman]?.  4], 
but  with  velocity  dispersion  as  in  the  theory  of 
Holden  and  Defines  [']. 

The  development  of  a  complete  noise  theory  for 
ballistic  and  near-ballistic  devices  is  being  planned: 
such  a  theory  may  aid  considerably  in  predicting  the 
mode  of  operation  (ballistic,  near-ballistic,  non- 
ballistici  from  the  measured  high-frequency  noise. 

ft.  DISCISSION  Of  n  '  pn  ’  RKStl.TS 

These  devices  showed  large  excess  noise.  The  noise 
is  not  very  close  to  I  /.  If.  nevertheless,  we  apply 
Hooge's  formula,  at  10  H/  and  100pA.  we  obtain 
h.  ^  .1  •  It)  ''. 

l  or  I  >  '00  mV.  the  I  I  characteristic  of  Fig.  5 
is  in  reasonable  agreement  with  the  theoretical  pre¬ 
dictions  and  previously  reported  results [1).  The  cur¬ 
rent  below  the  punch-through  voltage  (  i  1 50  mV)  is 
not  well  understood  It  is  significant  to  note  that  the 
character  of  the  noise,  in  particular  its  spectral  shape, 
does  not  change  when  we  pass  the  punch-through 
voltage,  see  fig  10.  only,  above  the  punch-thr ":gh 
volt.ire  the  noise  magnitude  starts  to  decline,  going 
no  loneer  as  /  .  see  f  ig  11  We  still  remark  that 
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similar-type  spectra,  going  slower  than  I  /.  were 
observed  in  bpm  p  'up  '  punch-through  diodes  by 
van  de  Roer[22|. 

van  der  Ziet[2?|  has  suggested  an  explanation  for 
the  n  '  pit  ’  d.c.  characteristics,  as  well  as  the  noise 
spectra.  He  interprets  the  low-frequency  noise  as  a 
form  of  Hooge  type  1  /  noise,  with  a  high  x„.  He 
notes  that  there  is  a  large  spillover  of  electrons  into 
the  p  region  where  they  recombine  with  the  holes. 
Though  there  are  relatively  few  holes,  they  control 
the  transport  properties  ambipolarly  Hence  there  are 
many  collisions  in  the  p  region  and  the  noise  is  high 
At  larger  forward  bias,  the  hole  recombination  in  the 
p  region  becomes  so  low  that  it  can  no  longer  control 
the  electron  transport  Thus,  the  number  of  collisions 
decreases  and  the  current  becomes  more  ballistic;  as 
a  consequence  the  (/.  I  )  characteristic  curves  upw  ard 
and  the  noise  decreases. 

The  validity  of  this  explanation  could  be  checked 
if  noise  measurements  were  made  on  p  '  np  '  and 
p  '  pp  *  devices.  In  the  p  ’  pp  ’  dev  ices  one  expects  a 
near-linear  (/.  L)  characteristic  and  high  noise;  this 
noise  should  rise  with  l\  and  not  go  through  a 
maximum.  The  />  ‘  np  '  devices  should  show  ambi- 
polar  transport  governed  by  electrons.  The  transport 
should  then  be  largely  ballistic  at  low  bias.  At  higher 
bias,  the  electrons  would  cease  to  act  as  ambipolar 
agents,  causing  the  transport  to  become  nonhallistic; 
the  (/.  !')  characteristic  would  go  slower  than  linear 
and  the  noise  u ould  rise  strongly,  i.e..  much  faster 
than  T .  Noise  measurements  on  p  '  pp  '  and  p  '  np' 
devices  arc  in  the  planning  stage. 

CONd.lSIONS 

Near-ballistic  n  '  nn  '  devices  exhibit  extremely 
low  1  /  noise,  with  a  Hooge  parameter  of  b.O  x  10  *. 
This  indicates  that  the  I  / noise  is  probably  caused  by 
lattice  scattering,  due  to  polar  optical  phonon  emis¬ 
sion.  which  is  rare  in  the  near-ballistic  regime.  The 
thermal  noise  of  these  devices  is  slightly  higher  than 
Nyquist  noise  for  the  highest  current  levels  observed. 
It  indicates  a  near-ballistic  origin,  affected  by  carrier 
drag  effects  and  by  the  velocity  dispersion  of  the 
injected  carriers. 

Near-ballistic  n  '  pn  '  dev  ices  exhibit  very  large 
low-frequency  noise  The  initial  current  and  the  noise 
may  be  caused  by  ambipolar  control  of  the  holes  in 
the  p  region.  The  noise  rises  initially  as  /:,  then  goes 
through  a  maximum,  and  then  decreases,  due  to  the 
ballistic  regime  overtaking  the  initial  current  regime. 
Thermal  noise  for  these  devices,  requiring  mea¬ 
surements  above  100  MHz.  have  not  yet  been  carried 
out. 

Similar  measurements  on  n  '  nn  4  diodes  arc  being 
reported  by  Peczalski  <  r  u/.  [24] 
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Since  the  appearance  of  the  above  article,  much  more  work  has  been  done 

■{•  •“ 

on  near-ballistic  devices,  in  particular  on  n  p  n  devices  (see  below)  and  on 

+  —  +  — 

n  n  n  devices  with  different  widths  for  the  n  layers  (see  the  next  section). 

^ 

The  n  p  n  device  characteristic  was  measured  at  different  temperatures, 
see  Fig.  1  of  Sec.  E.  There  is  sti  11  a  linear  1-V  region  at  low  bias,  then  there  is 
a  superlinear  region,  while  at  high  bias  the  different  temperature  curves  con¬ 
verge.  We  believe  that  the  latter  asymptote  represents  the  true  ballistic 
range.  The  results  appear  even  clearer  if  we  represent  the  data  in  a  plot  of 
DC  conductance  I/V  versus  the  bias  voltage  V.  This  is  shown  in  Figure  2.  We 
now  see  clearly  that  there  is  an  initial  low  conductance,  which  depends  on  the 
temperature  (though  below  77 K  the  data  seem  to  coincide);  then  there  is  a 
transition  region,  followed  by  a  final  temperature-independent  ballistic  range. 
The  noise  spectra  could  initially  only  be  measured  in  the  low  conductance 
region.  The  results  for  300 K  are  shown  in  Fig.  3  and  the  results  for  77 K  are 
shown  in  Fig.  4.  We  note  that  the  slopes  change  with  temperature.  The  300 K 
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spectra  are  of  the  form  1/f^'^  -  1/f^*^,  while  the  77 K  spectra  are  close  to  the 


form  1/f. 


The  low  bias  current  has  puzzled  many  investigators,  see  the  original 
papers  by  Eastman  et  al.  (e.g.,  M.S.  Shur  and  L.F.  Eastman,  Electronic  Letters 
16,  522,  1980).  The  expected  characteristic  is  as  in  the  dotted  line  in  Fig.  1 
The  bias  of  100  mV  (20  mV  per  layer  since  the  mesa  contains  five  devices  in 
series)  is  the  "punch-through"  voltage.  The  current  for  lower  biases — which 


should  not  be  there  according  to  the  pure  ballistic  theory — is  called  the  pre¬ 


punch  through  current.  Supported  by  the  noise  data,  we  believe  that  this 


current  is  recombination  current,  not  unlike  that  in  p-n  junction  depletion 
layers,  as  in  the  standard  theory  of  Noyce,  Shockley,  and  Sah.  First  of  all, 
we  note  that  for  zero  bias  and  for  very  low  bias  the  p  region  is  depleted  of 
holes,  due  to  spill-over  of  electrons  from  the  n+  regions,  as  was  calculated 
quantitatively  by  van  der  Ziel  and  Shur.  Thus  the  acceptor  sites  are  all  empty 
the  electrons  being  transferred  to  the  valence  band.  At  sufficiently  forward 


bias,  the  electrons  that  are  injected  will  be  captured  by  the  acceptors,  and 


eventually  by  the  valence  band.  Consequently,  there  is  recombination  current 
flowing,  until  all  acceptors  are  filled  (those  before  and  beyond  the  potential 
minimum).  This  process  is  apparently  accompanied  by  a  very  large  amount  of 
recombination  noise,  with  distributed  time  constants,  resulting  in  1/f  ,  a<l, 


noise.  For  low  temperatures  a  increases  since  the  shorter  time  constants  then 
become  less  prevalent.  This  noise  is  therefore,  in  our  opinion,  not  "true 
1/f  noise"  but  a  form  of  g-r  noise.  This  is  also  born  out  by  the  fact  that 
the  DC  conductance  in  the  prebreakdown  region  is  temperature  dependent.  Above 
the  "punch-through"  voltage  (all  acceptors  now  occupied,  n>N  )  the  ballistic 
current  appears.  The  DC  conductance  limit  of  Fig.  2  represents  a  conductance 
which  is  approximately  a  factor  three  lower  than  in  the  n+n"n+  samples.  This 


is  due  to  the  deeper  space  charge  minimum  of  the  former  device.  Note  that  the 
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space  charge  is  q(n-ND)  in  the  n-layer  while  it  is  (n  +  N^)  in  the  p-  layer. 
Bosnian  has  made  quantitative  computer  calculations,  based  on  the  Fry-Langmuir 
theory,  similar  to  those  by  Holden  and  Debney  (preceding  paper.  Ref.  7).  He 
found  the  computed  DC  ballistic  conductance  limit  to  be  in  excellent  agreement 
with  the  measurements,  both  for  the  n+n  n+  device  and  for  the  n+p  n+  device. 

A  paper  on  these  data  and  computations  is  in  preparation. 

Ill .  Further  measurements  in  GaAs  structures 

a)  0.24y,  0,4y,  and  l.lu  devices.  Jean  Andrian 

1.  Low-frequency  measurements  and  I-V  characteristic  of 
n+n  n+  diodes  with  different  lengths 
We  have  made  systematic  noise  measurements  on  n+n  n+  GaAs  structures 
for  various  temperatures  (77  K  and  300  K)  and  various  bias  currents.  These 
measurements  are  shown  in  Figs.  5,  6,  7,  8  and  9.  In  order  to  compare  the 
magnitudes,  we  also  plotted  (S-j-Zl2)*  f,  which  corresponds  to  a/N  in  Hooge’s  model 
of  1/f  noise,  as  a  function  of  the  bias  current.  This  is  shown  in  Figs.  10,  11 
and  12.  Next  we  plotted  this  quantity  as  a  function  of  the  length's  device. 

The  results  are  shown  in  Fig.  13.  As  we  note,  the  noise  of  the  0.27y  device  is 
higher  than  that  of  the  0.4y  device.  This  is  presently  not  understood,  since 
the  0.27y  device  is  the  nearest  to  collisionless,  ballistic  behavior. 

Besides  noise  measurement,  we  have  measured  the  I-V  characteristics  of 
these  devices.  These  are  shown  in  Figs.  14  and  15. 

The  0  •  4y  device  is  linear,  as  was  observed  by  Schmidt  et  al. ;  see  pre¬ 
ceding  section.  In  the  1.1  um  we  observed  a  change  of  conductance  at  higher 
bias  («0.7  volt).  This  phenomenon  does  not  occur  in  the  .4  urn  device  for  the 
same  applied  field.  We  believe  that  the  .4  ;.m  device  operates  in  the  near- 
ballistic  mode,  while  the  1.1  um  one  is  a  bulk  device.  This  change  of  con¬ 
ductance  is  probably  due  to  intervalley  transfer,  i.e.,  the  electrons  go  from 
a  valley  with  low  effective  mass  to  a  valley  with  a  higher  effective  mass. 


We  should  be  able  to  detect  the  noise  due  to  this  intervalley  scatter¬ 
ing.  This  would  be  a  very  interesting  phenomenon.  In  order  to  measure  that 
noise  we  set  up  a  noise  measurement  system  which  operates  in  the  radio  fre¬ 
quency  range.  With  this  system  we  will  be  able  to  measure  the  change  of  thermal 
noise  with  bias  and  the  intervalley  scattering  noise  in  the  transition  region. 

2 .  R-F  noise  measurement  system 

Essentially,  we  have  to  measure  the  thermal  noise  of  a  device  with  a 
lft  resistance.  We  set  up  the  correlation  measurement  system  indicated  in 
Fig.  16.  We  made  a  preliminary  test  on  the  linearity  of  the  mixer  and  the 
sensitivity  of  the  system.  We  plotted  the  DC  meter  reading  versus  the  power 
input  at  the  power  splitter.  The  results  were  excellent.  The  mixer  has  a  wide 
range  in  which  it  operates  linearly;  however,  we  need  to  have  a  DC  amplifier 
with  a  gain  of  1000  (linear)  in  order  to  be  able  to  observe  the  thermal  noise 
with  a  bandwidth  of  1  MHz.  This  is  done  by  a  single  invertor  circuit  as  shown 
in  Fig.  17.  Measurements  with  this  high-frequency  noise  system,  in  an  attempt 
to  find  the  intervalley  noise,  are  in  progress. 

-f.  _  -f- 

b)  Pulsed  noise  measurements  on  n  p  n  GaAs  devices.  Chris  Whiteside 

1.  Pulsed  bias  noise  measurement  system 

Completion  of  the  pulsed  bias  noise  measurement  system  has  been  accom¬ 
plished.  This  system  permits  noise  measurements  on  devices  at  very  high  bias 
levels.  The  device  is  pulsed  at  a  low-duty  cycle  in  order  to  avoid  excessive 
Joule  heating.  Figure  18  shows  the  final  pulsed  noise  measurement  system. 

2.  Results 

Initial  experiments  using  the  pulsed  system  were  completed  on  a  near- 
ballistic  GaAs  diode.  This  diode  had  been  used  in  previous  experiments.  We 
wanted  to  extend  the  bias  current  level  in  order  to  determine  the  current  noise 
in  the  ballistic  range.  According  to  the  I-V  characteristic  of  Fig.  1  (see 
previous  section),  it  is  believed  that  the  device  becomes  ballistic  at  high 
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current  levels.  Thus  a  drop  in  the  current  noise  must  take  place.  Fig.  19 

shows  that  as  the  current  level  is  increased,  the  noise  begins  to  fall  off. 

This  is  a  plot  of  vs.  f.  The  current  through  the  device  was  increased  up 

to  one  ampere.  Device  heating  was  kept  low  by  the  low-duty  cycle  of  the 

measurement  system.  Clear  plots  of  the  noise  as  a  function  of  current  for 

4  MHz  and  8  MHz  are  given  in  Figs.  20  and  21.  Though  the  decrease  of  the  noise 

is  very  pronounced,  we  have  not  been  able  to  reach  the  ballistic  noise  limit, 

-14 

which  is  estimated  to  be  at  10  at  the  ordinate  scale  of  Figs.  20  and  21. 

This  would  require  current  in  excess  of  10  amperes!  Measurements  as  a  function 
of  temperature  are  in  the  planning  stage. 

c)  Superlattice  measurement 

We  have  obtained  a  three-layer,  GaAs-Al  Ca,  As  heterostructure  device 

x  1— x 

(from  Dr.  Morkoc)  with  which  noise  measurements  can  be  made.  Heterostructures 

make  it  possible  to  measure  the  effects  of  a  two-dimensional  electron  gas. 

The  structure  of  these  devices  is  shown  in  Figure  22,  with  the  physical 

dimensions  of  the  layers  given  in  Table  1.  The  A1  mole  fraction  of  20%  causes 

a  conduction  band  discontinuity  of  0.21  eV.  This  bandgap  difference  between 

the  GaAs  and  A1  Ga,  As  layers  forms  potential  wells.  Electrons  confined  in 
x  1— x 

these  potential  wells  behave  like  a  two-dimensional  electron  gas. 

The  current-voltage  characteristic  of  the  heterostructure  is  shown  in 
Figure  23.  Noise  measurements  were  also  performed  and  are  shown  in  Figure  24. 
All  measurements  were  performed  at  300  K.  Problems  involved  with  cooling  the 
heterostructure  have  been  solved  so  that  noise  spectra  as  a  function  of  tem¬ 
perature  will  be  measured  in  the  near  future. 


Table  1 


Device  (three-period) 

#337 

#328 

GaAs  top  layer  (undoped) 

200  A 

200  A 

A1  Ga.  As 
x  1— X 

(Si-doped) 

250  A 

— 1018/cm^ 
x  =  0.20 

1000  A 

,„18,  3 

10  /cm 

x  =  0.20 

Separation  layer 
(undoped  Al^Ga^  ^As) 

150  A 

150  A 

_? 

Sheet  electron  concentration  ng[cn  “]  300  K. 

77  K 

10  K 

l'* 

2.8  101*- 
12 

2.24  10 

1 2 

2.24  lO1 

1? 

2.67  10 

12 

1.02  10 

12 

1.88  10 

tJHall^Cm'”//^  3ec^  300  K 

77  K 

10  K 

6.65  103 

9.3  104 
2.06  105 

6.09  103 

7.63  104 

1.64  105 

Contact  spacing  L 

400  um 

400  ym 

Contact  width  W 

40  pm 

40  pm 
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Abstract 

Counting  statistics  of  alpha  particles  from  Americium-241  were  deter 
mined  over  periods  from  1  minute  to  1,000  minutes.  In  particular,  the 
two-sample  variance  or  Allan  variance  was  determined  for  many  sample  runs. 
According  to  a  recent  theorem,  there  is  a  unique  relation  between  the 
particle  flux  spectral  noise  density  and  the  Allan  variance.  It  was  found 
that  for  small  counting  periods  the  statistics  were  Poissonian,  correspond 
ing  to  shot  noise  of  the  particle  flux.  For  long  periods  the  counting 
statistics  were  found  to  be  non-Poissonian,  and  a  flicker  floor  of  ®  10 
was  established.  Good  agreement  with  the  quantum  theory  of  1/f  noise  was 
obtained.  These  experiments  are  the  first  quantitative  confirmation  of 
this  theory. 
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1.  Introduction 

Whereas  most  of  noise  phenomena,  like  shot  noise,  thermal  noise,  and 
generation-recombination  noise,  are  well  understood,  1/f  noise  remains  to 
some  extent  an  enigma.  This  noise  has  been  observed  in  a  great  variety  of 
systems:  semiconductor  devices,  music,  traffic  flow,  hourglass  flow,  the 
frequency  of  sunspots,  the  light  output  of  quasars,  etc.  Because  of  its 
universality,  some  investigators  believe  that  there  must  be  some  common 
phenomenon  operative  in  all  of  these  manifestations^^ 

One  of  the  general  theories  of  1/f  noise  is  the  quantum  theory,  based 
on  infrared  divergent  coupling  of  the  system  to  the  electromagnetic  field  or 
other  elementary  excitations.  It  was  mainly  developed  by  Handel^^^^.  This 
theory  is  fundamental  in  the  sense  that  it  derives  the  1/f  spectrum  from 
basic  quantum  physics  at  the  level  of  a  single  charged  particle  subject  to 
scattering  with  small  energetic  losses  due  to  bremsstrahlung,  although  the 
final  result  depends  essentially  on  the  presence  of  many  carriers;  in  this 
respect  1/f  noise  is  similar  to  electron  diffraction  which  is  a  one-particle 
effect,  but  which  can  be  seen  only  if  many  particles  are  diffracted.  In 
addition,  the  theory  is  universal  in  the  sense  that  any  infraquanta  with 
infrared-divergent  coupling  to  the  current  carriers  will  give  a  contribution 
to  the  observed  1/f  noise  proportional  to  their  coupling  constant.  Such 
infraquanta  are,  for  example,  very  low-frequency  photons,  various  types  of 
phonons,  shallow  electron-hole  pairs  on  the  Fermi  surface  of  a  metal,  spin 
waves,  correlated  states,  etc. 

Recently,  this  theory  has  been  reformulated  with  quantum-optical 
terminology  and  compound-Poisson  statistics  in  a  paper^^  written  by 


Van  Vlict  and  Handel,  which  led  to  the  idea  of  verifying  the  theory  on  a 
"clean"  system  outside  the  domain  of  currents  in  solids:  radioactive  a-decay. 


In  this  paper  we  present  data  obtained  from  extensive  measurements 
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on  counting  techniques  for  a-particles  radioactive  decay  from  Am,  which 
have  shown  that  the  statistics  are  non-Poissonian  for  large  counting  times 
(order  1,000  minutes)  in  contrast  with  the  fact  that  many  textbooks  cite 
a-decay  as  the  example  "par  excellence"  for  Poisson  statistics.  With  the  help 
of  the  newly  devised  Allan  variance  transform  theorem^ we  found  that  this 
excess  noise  does  have  a  1/f  spectrum,  and  the  "flicker  floor"  due  to  the 
presence  of  1/f  noise  in  the  decay  has  a  value  of  10  ^  which  is  in  fair  agree¬ 
ment  with  Handel's  quantum  1/f  noise  theory.  This  result  indicates  that  1/f 
noise  is  caused  by  emission  of  long  wavelength  infraquanta,  such  as  soft 
photons  causing  minute  inelastic  losses  in  the  scattered  wave  packet. 

2 .  Results  of  the  Theory 

2.1.  Quantum  approach  to  1/f  noise 

It  is  known  that  upon  scattering  a  beam  of  electrons  will  emit 
bremsstrahlung.  The  power  spectrum  W(f)  of  the  emitted  radiation  is  indepen¬ 
dent  of  frequency  (W  =  constant)  at  low  frequencies  and  decreases  to  zero  at 

an  upper  frequency  limit  f  which  is  approximated  by  E/h,  where  E  is  the 

m 

kinetic  energy  of  the  electrons;  h  is  Planck's  constant.  Consequently,  the  rate 

W 

of  photon  emission  per  unit  frequency  interval  is  N(f)  =  i.e.,  propor¬ 

tional  to  1/f.  Therefore,  we  conclude  that  the  fraction  of  electrons 
scattered  with  energy  loss  c  is  proportional  to  1/e,  i.e.,  the  relative 
squared  matrix  element  for  scattering  with  energy  loss  e  is  |b^,(e)|  cxl/e. 

If  the  incoming  beam  of  electrons  is  described  by  a  wave  function 
exp  [  (i/ii)(p*r-  Et)  ]  ,  the  scattered  beam  will  contain  a  large  nonbremss t rahlung 
part  of  amplitude  a,  and  an  incoherent  mixture  of  waves  of  amplitude  ab^(c) 

with  bremsstrahlung  energy  loss  c  ranging  from  some  resolution  threshold  c  to 

o 

an  upper  limit  A<E,  of  the  order  of  the  kinetic  energy  E  of  the  electrons 


r  _  ,  /  A 

i|jt  =  exp  J^l^jCpT  -  Et)  a  1  +  J 

e 


bT(e)ei£t/^  d A  . 


(2.1) 


Here  b^,(e)  =  |b^,(e)|e  has  a  random  phase  y  which  implies  incoherence  of 
all  bremsstrahlung  pares,  and  Jb  (e)]‘"  is  proportional  to  1/c,  as  we  saw 
above. 

In  Eq.  (2.1)  the  f requency-shif ted  components  present  in  the  integral 
interfere  with  the  elastic  term,  yielding  beats  of  frequency  e/fi.  The 
particle  density  given  by  Eq.  (2.1)  is 


i«Tr  =  I a 


2l(1  +  2  I  COS  (^r  +  Ye)  dc 


A  n 

i  /■■' 


(c)bT(c’)el(e''e)l:/ri  dede'  ; 


o  O 


(2.2) 


the  second  term  in  large  parentheses  describes  the  particle  density  beats. 

,  ,2  ,  ,  2  ,  <  2 


autocorrelation  function  will  be 


~  2  a 


b  ( c )  1  cos 


i  is  defined  by  5 | 

>|2  =  1 

M2  -  < 

t  1  ^  t+x  >  ~  ^  1 

2  2 

> 

(¥)dc  =  2  ]  a  |  4  ^ 

A/h 

/  * 

|b(c) 

f 

o 

o  o  (2.3) 
which  is  proportional  to  jb(c)|“  and  hence  proportional  to  1/f.  Therefore, 
the  spectral  density  of  the  particle  concentration  fluctuation  (the  Fourier 
transform  of  Eq .  (2.3)  )  is  proportional  to  1/f. 

The  relative  bremsstrahlung  rate  ;b(;)j‘  can  be  derived  as  follows. 


The  constant  spectral  energy  density  can  be  written  as  w(f)  =  juvj  ”V3c-'<, 


where  e  is  electronic  charge,  *:  the  dielectric  constant  of  the  medium. 


r 
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c  is  the  velocity  of  light,  and  Av  is  the  velocity  change  in  the  scattering 
process.  The  relative  scattering  rate  density  with  energy  loss  z,  jb(e)|2, 
is  obtained  by  dividing  W(f)  by  the  energy  of  a  photon  e  =  hf. 


where 


|  b  ( f )  |  2  =  =  |£  liMl 

3c3<2irfif  *C  37IKf 

|b(e)|2  =  =  oa 

3c3Kh£ 


.  _  2(A8)"  Av 

A  =  - -  *  AS  =  —  , 


(2.4) 


(2.5) 


(2.6) 


e2  e2 

and  a  = jp—  is  the  fine  structure  constant.  In  the  M.K.S.  system  a  =  - 

ui c  J  die  < 

o 

where  is  the  dielectric  constant  of  vacuum.  The  spectral  density  of 
the  relative  fluctuations  is  from  (2.3),  (2.5),  and  the  Wiener-Khintchine 
theorem, 

S  |  1 2(f)  /<  j  :-|V  =  211+  :uUn(A/c  )f2aA/<f  *  2aA/<f  .  (2.7) 

|  v  |  o 

2.2.  The  Allan  variance  transform  theorem 

The  main  link  between  counting  statistics  and  particle  current  noise 

r  g  I 

is  provided  by  MacDonald's  theorem1  , 


CO 

<  AM-p  >  =  J  sm(.J^  1sintoT  d 


(2.S) 


witli  inversion 


S  (.0 
m 


CO 


(2.9) 


here  <  ,'.M^  >  is  the  variance  of  the  total  number  of  particles  detected  in 

a  time  interval  (t,  t + T)  and  S  (  )  the  noise  spectral  density  of  the  flux 

m 


r 


fluctuations  'm(t).  This  theorem  is  useful  for  Poissonian  statistics. 
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Unfortunately ,  for  1/f  noise  Eq.  (2.8)  is  not  applicable,  since  the  integral 

diverges.  However,  a  useful  concept  in  this  case  is  the  "two  sample  variance" 

or  "Allan  variance" ^ ^ ^  ,  Let  he  the  average  counting  rate  in  (t,  t  +  T) 

and  m^,2^  the  counting  rate  in  (t  +  T,  t+2T).  Likewise,  let  be 

(2)  (9) 

the  total  number  counted  in  (t,  t  +  T)  and  let  =  m2  T  be  the  total  number 

counted  in  (t  +  T,  t  +  2T) .  Then  the  Allan  variance  is  defined  by 


A2 

“T 


< 


(D  (2) 
mT  ~mT  I  > 


(2.10) 


and 


>=rV; 


(2.11) 


A2  A  -> 

The  variance  a  “  (which  means  (a  )~)  turns  out  to  be  finite  for  1/f  noise. 


The  transform  theorem  reads 

“  S(c) 


[7] 


■  T 


I 


■  4  |  -)T 

cm  I  - 


sin 


d.j 


(2.12) 


with  inversion 


S  (w)  =  - 
m 


i"+£ 


AJ 


dn 


1 

cos  pn 


uT2  1  -  2P_3 


r(p) 


f  dT  A2 

J  TP  \ 


(T) 


(2.13) 


-  i>»-f  8 

(Slightly  more  complicated  inversion  forms  using  partial  Mellin  transforms 

are  found  in  [7].)  For  Poissonian  shot  noise  S  (w>)  =  2m  ,  where  ra  is  the 

mo  o 

average  counting  rate.  Substituting  S^Cj)  into  Eq .  (2.12),  one  has 

A 2"C 
c  “(T)  =  m  T.  For  1/f  noise,  with  a  spectrum  of  S  ( .)  =  r,  where  C  is  a 
o  m  ,  | 

1  A9 

constant,  the  Allan  variance  J  =  2CT"in2.  Now  suppose  that  the  noise  is 

T 

composed  of  shot  noise  and  1/f  noise,  i.e.. 


aJ2(T)  =  m  T  +  2CT‘  :.n2  ; 


(2.14) 
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We  recall  that 
variance  R(T)  = 


<Mr> 


m  T,  so  that  a  measurement  of  the  relative  Allan 
o 


4^(T)  /<Mt>2  yields, 


R(T)  =  -i=  +  2C’in2  ,  (2.15) 

m  i 
o 

C  * 

where  C'  =  — —  is  the  characteristic  strength  of  the  1/f  noise  S  (f)/m“  . 

m  -  m  o 

For  short-time  intervals  the  term  — is  dominant,  hence  R(T)  is  proportional 

‘  o 

to  1/T.  When  T  is  long  enough,  the  term  2C'i,n2  becomes  dominant.  For  large  T 
one  cannot  further  reduce  the  relative  accuracy  by  longer  counting; 

R(T  )  =  2C'in2  is  therefore  called  the  "flicker  floor". 
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We  conclude  from  the  above  that  the  presence  of  1/f  noise  in  counting 
statistics  can  be  determined  from  a  measurement  of  the  Allan  variance  as  a 
function  of  T. 


3.  Experimental  Method 
3.1.  Procedure 

The  block  diagram  of  the  counting  system  being  used  to  investigate  1/f 

fluctuation  in  the  a-particle  emission  rate  is  shown  in  Fig.  25.  The  source 
.  241 

is  An,  which  decays  with  a  half-life  of  T^7  =  458  years  with  the  emission 
of  5.48  (86%),  5.44  (12.7%),  and  5.34  (1.3%)  MeV  a-particles  into  “^Np.  The 
detector,  a  silicon  surface-barrier  detector,  is  reverse  biased  at  80  volts, 
and  the  dead  times  of  the  ND575  Analog  to  Digital  Convertor  (ADC)  and  ND66 
Multi-Channel  Analyzer  (MCA)  are  60  nsecs  and  6  usees  respectively.  There¬ 
fore,  no  dead-time  correction  is  necessary,  as  long  as  the  counting  rate  is 

f  10 1 

kept  below  1,000  counts  per  second  (or  the  averaged  time  elapse  between 

two  counts  is  higher  than  1,000  usees). 

A  typical  full  energy  spectrum  measured  in  these  experiments  is  given 
in  Fig.  2b,  in  semi- 1 ogar i thnie  scale.  The  spectrum  is  shown  on  a  displav 
screen  while  accumulating  counts  and  the  final  result,  after  a  chosen  time  T, 


is  stored  in  the  histogram  memory  units  of  the  ND66  MCA.  The  full  width 


half  maximum  (fwhm)  of  the  spectrum  can  be  found  by  moving  the  CRT  display 
cursor  across  the  alpha-peak  channel  regions,  with  the  cursor  display  indi¬ 
cating  the  number  of  counts  within  each  energy  interval  (channel) .  The  sys¬ 
tematic  range  used  in  these  experiments  with  respect  to  the  5.48  Mev  peak 
channel  was  from  six  fwhm's  below  the  peak  channel  to  two  fwhm's  above  the 
peak,  in  this  manner  spanning  the  three  prince  peak  radiation  types.  There¬ 
fore,  with  the  range  being  a  function  of  the  detector  resolution,  the  total 
number  of  counts  M^,,  which  will  be  analyzed  later,  is  always  a  fixed  portion 
of  the  full  spectrum. 

The  counts  M^,  of  adjacent  time  intervals  can  be  read  directly  from  the 
memory  units  of  the  ND66  MCA;  thus  the  Allan  variance  can  be  calculated  by 


(3.1) 


since 


<\>  -T  X,  >l 

1=1 


(i) 

T 


(3.2) 


the  relative  Allan  variance,  R(T) ,  defined  by  Eq.  (2.15),  can  be  found  by 
using  Eqs.  (3.1)  and  (3.2). 

Part  of  the  data,  mostly  the  total  number  of  counts  for  T  longer  than 


1000  minutes,  were  not  read  directly  from  the  ND66  Multichannel  Analyzer. 


.(1) 


An  "add-up"  method  was  adopted;  namely,  M^OOO  was  Stained  by  adding  up 


M-i.  and  M-qq.  and  M^qqq  equals  the  sum  of  and  M-qq>  etc.  Physically, 


500 


1000 


since  M^qqS  were  measured  in  adjacent  time  intervals,  of  course  the  first  two 


can  be  added  up  as  the  total  number  of  counts  for  the  first  1000-minute 


(2) 
Mv  ' 


interval,  and  the  third  and  fourth  can  be  summed  up  as  M^qoo'  Howcver,  in 


order  to  check  the  validity  of  this  method,  the  following  experiment  has  been 
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By  making  use  of  a  T-connector,  the  output  signal  of  the  Ortec  410  ampli¬ 
fier  was  fed  simultaneously  into  two  ND575  ADCs  (see  Fig.  25).  The  first  one 
(ADC  ill)  counted  100-minute  measurements  for  310  times,  and  the  second  one 
(ADC  #2)  accumulated  500-minute  counts  for  62  times;  hence,  both  ADCs 
covered  exactly  the  same  time  span.  The  "add-up"  method  was  applied  to 

M,„„s,  obtained  from  ADC  #1,  to  find  out  the  calculated  M_r  s . 
iUU  ouu 

Table  1  lists  the  results  calculated  from  both  ADCs  for  500-minute 
measurements.  The  difference  between  them,  in  each  category,  is  <  1%.  This 
shows  the  validity  of  the  "add-up"  method. 


ADC  ill 


ADC  it 2 


<M500>  9062A6A.  36 

</2  .  7862171.3 

H500 

R(500)  9.573  10~ 


9067009.84 
7943285.5 
9.662  10' 


Table  1.  Comparison  of  the  results 
obtained  from  the  "add-up”  method  (ADC  #1) 
and  the  real-time  measurements  (ADC  #2) . 


Because  of  the  existence  of  "variance  of  variance"  (or  "variance  noise") , 
the  relative  Allan  variance  itself  is  a  fluctuating  parameter.  In  order  to 
obtain  an  accurate  value  of  R(T) ,  a  sufficient  number  of  measurements  must  be 
made,  especially  when  T  is  short.  Figure  27  shows  R(T)  versus  the  number  of 
measurements,  N,  for  T  =  1  rain.  When  N  is  small,  R(T)  is  spread  over  a  wide 
range.  When  N  is  increased,  R(T)  shows  less  spread  and  finally  converges  to 
a  stable  value.  The  minimum  value  for  N’  is  about  70.  Figure  28  gives  a 
similar  plot  for  T  =  3  rain.  For  a  reliable  value  of  R(T) ,  N  has  to  be  larger 
than  50. 

Due  to  the  presence  of  variance  noise,  the  experimental  results  of  R(T) 
contain  a  fluctuation  term  .R,  i.e.. 
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(3.4) 


R  (T)  =  — -  +  2C'Jln2  +  AR  ,  (3.3) 

exp  m  T 

o 

where  the  average  value  <  AR>  should  be  zero.  The  experiment  was  then 

repeated  several  times,  and  the  average  of  R  (T) , 

exp 

<  R  „  (T)  >  =  +  2C'in2  +  <AR>  ,  (3.4) 

exp  m  T 

o 

is  obtained.  Although  the  value  of  <AR>  in  Eq.  (3.4)  can  hardly  be 

exactly  zero  for  a  finite  number  of  measurements,  it  should  be  reduced  a 

great  deal  compared  with  the  value  of  AR  for  single  measurement.  The  value 

of  <R  (T)  >  gives  the  best  estimation  to  the  true  value  of  R(T) . 
exp 

The  distance  between  radioactive  source  and  detector  was  adjusted 
such  that  very  close  counting  rates  were  obtained  while  repeating  measurements. 
However,  it  is  difficult  to  obtain  identical  counting  rates.  Therefore,  the 
shot  noise  term  in  relative  Allan  variance,  R(T) ,  is  slightly  different  betweeen 
each  series  of  measurement.  * 

Before  the  average  of  R  (T)  is  taken,  the  shot  noise  term,  -^r  , 

exp  m  T 

o 

should  be  normalized  to  the  same  counting  rate.  Here  a  rate  of  18,000  counts 

per  minute  was  chosen.  The  normalized  relative  Allan  variance  R  (T)  is  then 

n 


Rn(T)  Rexp(T)  '  ra  T  +  18,000  T  18,000  T  +  2C  Zn2‘ 
o 


(3.5) 


The  value  <Rn(T)>  is  now  used  to  estimate  the  true  value  of  R(T) .  This 
value  is  calculated  as  follows: 


<  R  (T)  >  = 

n 


£R  . (T)  x  DF. 
ni  i 


(3.6) 


where  R  . (T)  is  the  value  of  R  (T)  obtained  from  the  i  series  and  DF . 

ni  n  i 

stands  for  the  degrees  of  freedom  of  that  particular  value,  which  equals 


r 


the  number  of  measured  time  intervals  minus  one. 
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3.2.  Results 

Fig.  29  shows  the  experimental  data  for  <R^(T)>.  The  theoretical 

curve  (solid  line)  is  based  on  R(T)  =  (1/18,000  * T)  +  1  *  10  '  versus  1/T, 

which  suggests  that  the  value  of  2C'2n2  is  about  l  x  10  ^ . 

Fig.  30  shows  the  average  value  of  the  normalized  Allan  variance, 

< >  ,  which  is  obtained  by 
ILn 


A 0  0 

<0*“  >  =  <R  (T)  >  X  (18,000  XT) *■ 
Mj.n  n' 


(3.7) 


-7 

The  theoretical  curve  (solid  line)  is  based  on  a  “  =  18,000  T  +  1  *  10  x 

2  “t 

(18,000  *  T)  .  Very  good  agreement  is  obtained  as  shown  in  these  figures. 

The  value  of  <  R  (T) >  and  1/18,000  *  T  (shot  noise  level)  are  listed 
n 

in  Table  2.  It  shows  clearly  that  at  this  counting  rate,  18,000  counts  per 
minute,  for  T  longer  than  100  minutes,  1/f  noise  becomes  noticeable.  For  T 
longer  than  1000  minutes,  1/f  noise  totally  dominates  the  noise  spectrum. 


T(min. ) 

1 

2  5 

10 

20 

50 

<  R  (T)  > 
n 

337.4 

261.0  123.3 

53.66 

27.86 

11.56 

1 

18,000  T 

355 . 6 

2  7.8  111.1 

55.56 

27.78 

11.11 

T  (min.) 

100 

200 

500 

1000 

2000 

3000 

<  R  (T)  > 
n 

7.136 

3.825 

1.799 

1 .  t>4S 

1.159 

1.816 

1 

18 ,000  T 

5.556 

2.778 

1.111 

0.5556 

0.2773 

0. 1S52 

Table  2.  The  values  of  <R  (T)  >  and  1/18 ,000  *  T; 
-  _  n 

all  values  x  10  ' . 
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4.  Discussion  and  Conclusions 

4.1.  Possible  systematic  sources  of  error 

Besides  the  fluctuation  of  emission  cross  section  described  in 
Handel's  theory,  two  other  possibilities  may  contribute  to  the  fluctuation 
in  the  total  number  of  counts.  First,  the  pulse  height,  which  is  produced 
when  an  a-particle  is  absorbed  by  the  detector,  can  be  a  strong  function  of 
the  bias  voltage  applied  to  the  detector^"^.  If  the  bias  voltage  were 
unstable,  changes  in  the  bias  voltage  could  result  in  changes  in  the  deple¬ 
tion  depth,  varying  the  sensitive  path  length  experienced  by  the  a-particle, 
and  the  depletion  layer  capacitance.  Thus,  a-particles  with  the  same  energy 
would  have  been  registered  under  different  energy  channels,  and  would  intro¬ 
duce  unwanted  fluctuations.  To  minimize  such  fluctuations,  the  detector 
control  unit  selected  for  use  in  these  experiments,  Ortec  Model  210,  has 
quite  good  stability:  bias  voltage  variation  with  the  line  voltage  is 
<  ±0.005%  for  105  -  125  V  AC  input,  and  the  stability  is  ±0.01%^“^.  Also,  the 
use  of  a  charge-sensitive  preamplifier  (Ortec  Model  109)  diminished  the 
system's  sensitivity  to  changes  on  detector  capacitance.  Secondly,  since 
a  fixed  portion  of  the  full  spectrum  is  counted,  whether  particles  with 
energies  near  boundaries  fall  within  or  without  the  region  they  bounded  con¬ 
tributes  fluctuations  to  the  total  number  of  counts.  However,  basing  the 
range  of  spectrum  integration  on  multiples  of  the  fwhm  reduces  variation  due 
to  peak  spreading  and  gain  shift. 

Therefore,  a  wider  energy  range  should  have  less  fluctuations  in  the 
total  number  of  counts  in  a  certain  time  period.  This  is  indeed  the  case, 
see  Table  3.  However,  for  T  =  100  minutes,  the  difference  between  R(T)s  for 
wider  range  and  narrower  ranee  is  negligible.  For  T  =  300  minutes,  the  dif¬ 
ference  is  4%,  but  the  value  of  <R  (r'00)>  is  62%  higher  than  the  shot  noise 
level  (see  Table  2).  For  T  =  S000  minutes,  the  difference  is  11%,  and  the 
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value  of  <R  (3000) >  is  ten  times  the  value  for  the  shot  noise  level.  Hence 
n 

the  energy  range  chosen  in  these  experiments  affected  very  little  the  final 
results . 


100  rain. 


Ch.  127-147  Ch.  112-157 
(Peak  Channel:  142,  FWHM:  2.5  channels) 


<M1qo>  1812492.94 


A2 

ctm 

100 

R(100) 


2430286.7 


7.398  x  io 


-7 


1828630.89 

2472416.52 
7.394  x  10-7 


500  min.  <  M_ > 


A  2 


>1 

500 
R(  500) 


9067009.8 


79432S5 . 5 


9.662  x  io 


-8 


9147489.02 

7762255.74 
9.277  x  10-8 


3000  min .  <  M^qq  > 


54398712.1  54881014.1 


3000 
R( 3000) 


698293362 


2.360  x  io 


-7 


634854514 


2.108  x  10 


-7 


Table  3.  Comparison  of  experimental  results  for 
narrower  energy  range  (peak  channel  -  FWHM  x  6  to 
peak  channel  +  FWHM  x  2)  and  wider  energy  range 
(peak  channel  -  FWHM  x  12  to  peak  channel  +  FWHM  x  6. 


4.2.  Comparison  with  the  theory  of  quantum  1/f  noise 

According  to  Handel's  theory,  the  constant  C*  in  the  1/f  noise  term  in 
the  relative  Allan  variance  is,  cf  Kq .  (2.7), 

C  =  8  uV./k  (4.D 


wiie  re  is 


a 


coiierence  factor;  for  a-particles  it  is  expected  to  be  close 


35 


to  one;  a  is  the  fine  structure  constant  1/137,  and  A  =  2(Av)-/37tc2  where 

-*  [5] 

Av  is  the  velocity  change  of  the  particles  in  the  emission  process  , 

c  is  the  velocity  of  light.  For  a-particles  one  finds 


—  =  33.28  x  lCf7  x  —  (4.2) 

K  K 

where  E  is  energy  in  MeV  and  k  is  the  dielectric  constant  of  the  radio¬ 
active  material. 

In  these  experiments  E  =  5.48  MeV,  so  that  the  value  of  2C'i,n2  (the 
value  of  the  flicker  floor,  F)  is 


F  =  44 


33.28  x  10 


5.48 


x  £n2  =  505.6  x  10  7  x— 


(4.3) 


To  the  authors'  knowledge,  nobody  has  every  measured  the  dielectric 
constant  of  Am09,  the  a-particle  source  used  in  these  experiments;  thus  the 
value  of  k  in  Eq.  (4.3)  is  unknown.  However,  since  the  atomic  structure 
of  AmO,,  is  similar  to  that  for  U0? ,  the  dielectric  constant  of  UO^ 

(20.4  ±1.5^^,  21.7  ±0.5^4^,  and  21.0  ±1^"^  can  then  be  used  as  a  refer¬ 


ence. 

If  one  considers  the  following  factors  that  1)  the  resulting  aA  in 
Eq.  (4.2)  is  the  sum  of  contributions  for  all  types  of  infraquanta  partici¬ 
pating  in  the  energy  transfer,  and  maybe  only  part  of  them  were  detected, 
i.e.,  the  actual  value  for  ocA  in  these  experiments  may  be  smaller  than  that 
given  in  Eq.  (4.2);  2)  the  coherence  factor,  4,  is  always  less  than  unity; 

3)  the  dielectric  constant  of  Am0o  may  be  larger  than  21:  then  the  measured 
value  of  F  =  1  <  10  7  is  in  the  right  ballpark  to  verify  Handel's  theory. 

In  particular,  for  4  =  0.05,  •-  =  25,  one  obtains  F  =  1.01  v  10  7 ,  in  accord 
with  the  observed  value  of  Figs.  29  and  30,  We  believe,  therefore,  that 
these  experiments  constitute  the  first  experimental  verification  of  Handel's 
quantum  L/f  noise  theory. 
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In  conclusion,  a-particle  decay  has  statistics  ttfhich  for  large  counting 

3 

intervals  (^10  min)  are  non-Poissonian.  We  established  the  presence  of  a 
flicker  floor  R(Tot)  =  1.0  *  10  In  the  frequency  domain,  this  indicates  the 

presence  of  1/f  noise  in  the  particle  flux  for  frequencies  of  the  order  of 

-A 

10  h2  or  less.  The  magnitude  of  the  flicker  floor  can  well  be  explained  as 
electromagnetic  quantum  1/f  noise,  providing  C  =  0.05  can  be  explained. 

(Later  note:  We  presently  think  the  factor  K  =  25  does  not  belong  in 
the  formula.  Thus,  in  reality,  the  theoretical  noise  is  a  factor  500  too  high. 

See  our  NB  at  the  end  of  Section  V  of  this  report  (page  51).) 
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V.  Partition  1/f  noise  in  pentodes.  C.J.  Hsieh  and  A.  van  der  Ziel 

It  is  well  known  that  pentodes  show  partition  noise,  in  addition  to  space- 
charge  suppressed  shot  noise.  It  is  less  known  that  a  pentode  also  shows  par¬ 
tition  1/f  noise.  The  effect  was  first  observed  by  Schwantes  and  van  der  Ziel, 
Physica  2b_,  1143  and  1157  (1960).  This  partition  noise  can  be  isolated  by  a 
feedback  procedure  that  reduces  the  regular  flicker  noise  with  a  large  resistor 

R  in  the  cathode.  The  then  remaining  noise  is  partition  1/f  noise.  This 
c 

noise  was  remeasured  in  an  attempt  to  compare  it  with  Handel's  theory, 
a)  Quantum  picture  of  flicker  noise 

According  to  Handel's  theory! 10],  there  is  a  quantum  approach  1/f  noise. 
There  should  be  a  quantum  1/f  noise  part  (Vi2)  present  in  partition  noise  in¬ 
volving  charged  current  carriers.  The  carriers  also  have  infrared  divergent 
coupling  to  infraquanta.  The  spectral  density  of  the  relative  fluctuations  is 


< (An) ->f  _  ?  «A 

f  ’ 

<n> 

The  infraquanta  we  are  going  to  consider  here  are  photons,  with 


(1) 


c*A  =  -1^ 
3tt 


V  —  V 

final  initial 


(2) 


where  a  =  1/137  is  the  fine  structure  constant  and  vf^na^  is  the  velocity  of 

a  carrier  at  the  end  of  the  process;  v.  .  .  .  is  the  velocity  before  the 

ini l iui 

process • 

In  the  case  of  a  pentode,  three  grids  are  present,  but  the  screen  grid 
is  the  only  one  carrying  a  positive  potential.  Then,  the  three  independent 
noise  sources  present  in  a  pentode  include  the  partition  noise  S^(f),  the  anode 


noise  .9^(0,  and  grid  noise  S^.  The  spectral  density  are  given  by 
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aA 

S  =  T22el  +  2  — ~  (r2I  )?  (3) 

3  3  L  3 

aA 

S  =  r22el  +  2  — (r2I  )2  (3a) 

g  g  f  g 

aA 

S  =  (1  -  r2) 2el  +  2  -7s-  [(1- r2)I  ]  .  (3b) 

P  g  f  g 

Here  we  have  denoted  by  F2  the  space  charge  noise  suppression  factor  and 

A  =  1.224  x  10-5  (u  -  u  )  (4) 

ca  a  c 

A  =  1.224  x  io-5  (u  -  u  )  (4a) 

eg  g  c 

A  =  1.224  x  io-5  u  (4b) 

g 


where  u  ,  u  ,  u  are  the  potentials  of  the  anode,  the  cathode,  and  the  screen 
a  c  g 

<•) 

grid,  while  u  gives  the  velocity  of  the  smoothed  electrons  (fraction  1-F"),  when 
they  arrive  at  the  screen  grid.  Note  that  in  order  to  take  into  account  in¬ 
coherence  effects  properly  we  have  to  replace  f  by  Nf  in  all  previous  equa¬ 
tions  in  the  denominator.  There  N  is  the  average  number  of  electrons  involved 
simultaneously  in  the  process  considered  (e.g.,  transition  from  the  cathode  to 


the  grid  in  the  nonsraoothed  wav  for  «A  /Nf) . 

eg 

From  the  current  generator 

defined  above,  we  can  easily  obtain  the  spectral  density  of  noises. 

ST  (f)  =  S  (f)  +  S  (f) 
a  P 

(5) 

sT  (f)  =  S  (f)  =  S  (f)  +  S  (f) 
tg  e  e  P 

(3a) 

ST  (f)  =  S  (f)  +  S  (f)  . 

1  d  g 

(5b) 

For  the  partition  noise  Handel  presents 

the  expression 

2-»A 

s  (f)  =  —f*  (!-■';’)  I- 

p  f  c  R 

(6) 
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where  A  =  1.224  x  10  u  and  u  is  the  velocity  of  the  smoothed  electrons 
g 

( fraction  1-F^).  For  this  expression  Dr.  van  der  Ziel  expects  instead 


S  (f) 
P 


2aA 

_ g. 

f 


I  I 
a  g 


(7) 


The  argument  is  as  follows.  Consider  two  pentodes,  one  with  a  screen  grid 

transparency  X  and  one  with  a  screen  grid  transparency  ( 1  —  A ) .  Then  I'  =  I9 

a  z 

and  I*  =  I  .  The  partition  noises  and  the  partition  flicker  noise  should  be 

z  <1 

the  same,  so  (f)  should  not  change  if  and  I  are  interchanged,  and 
P  8 

Equation  (7)  accomplishes  this.  Moreover,  all  electrons  arriving  at  the 
screen  grid  should  distribute  and  not  only  the  fraction  (l-T*")*  One  might  also 
question  the  factor  A^;  A^  is  proportional  to  the  anode  potential  u^,  that  is 
to  the  energy  with  which  the  electrons  arrive  at  the  anode;  Acg  is  proportional 
to  the  screen  grid  potential  u^.  So  van  der  Ziel  thinks  it  would  be  better  to 
write 


Sp(£)  = 


lyJA  A 
»  ca  cs 


I  I  . 

a  g 


(8) 


In  that  case, 


2a-v/A  A 

s r  T  (f)  - - (r4  - 1)  I  I 

Ia,!g  f  c  a  g 


(9) 


so  that  the  fluctuations  in  I  and  I  are 


a  g 

saturated  (f^  =  1).  This  looks  intuitive 


fully  uncorrelated  if  the  cathode  is 
ly  right. 


b)  Feedback  effect  of  the  cathode  resistor 

We  apply  feedback  in  the  cathode  lead  by  inserting  R  ,  as  shown  in 
Figure  31.  Then, 


v 

g 


-  (i  + 


rat  g 


v  )  R 

K  C 


(10) 
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Simplifying  the  above  equation,  one  obtains 


-i  R 

_ c  c 

Vg  =  1  +  gmtRc 


so  that  the  noise  current  of  the  cathode  with  feedback  is 


i'  =  i  +  g  v 
c  c  mt  g 


=  i  LJai.1  =— i 

S  C  l  1  +  8mtRc  ]  1  + 


8mtRc 


Now  we  assume  that  the  ac  signals  distribute  between  screen  grid  and  anode 
in  the  same  way  as  i^.  In  this  case 

gma  =  gmtX  ;  gm2  =  gmt(1'A)  *  < 


With  feedback  resistor  R  ,  then. 


i*  =  i  + 
a  a 


gv=Xi+i-g  — 
ma  g  c  p  bma  1 


i  R 

c  c  ,  .  t  ,  . 

— - =  Xi'  +  a 

+  g  R  c  p 

mt  c 


So  the  spectrum  of  the  anode  with  feedback  is 


Xa  =  *c 


1p  ’1c=l+gR 
mt  c 


SI  V 

S  (f)  =  - 2 - -  +  S  (f) 

a  (1  +  g  R  )‘  p 

mt  c 


(15a) 


and  we  can  also  write  Equation  (15)  and  (15a)  in  terms  of  equivalent  current. 


i'-  =  2qAf  - - - r  +  I 

3  L(1  +  gmtRc>^  V 


(15b) 


I  =  - - - -  +  B 

(1  +  gmCRc)' 


(15c) 


where  A  and  B  arc  constants.  We  now  define  the  new  pentode  transconductance 


with  feedback 


41 


mt 


3mt  1  +  g  R 
mt  c 


(16) 


gmrX 

mt 


&ma  1  +  g  R  * 
mt  c 


(16a) 


i  does  not  change  since  it  flows  from  screen  grid  to  anode.  From  Equa¬ 
tions  (12),  (14),  and  (15)  we  can  get 


g2  X 

4KTR'  df  g  2  =  4KTR  df  - — 

na  ma  nc 


(1  +  g  R  ) 
mt  c 


+  4K.TR  df  g  “ 
np  ma 


(17) 


since 


2  \2 


.2  _  mt 


ma 


(1  +  VV 


Dividing  Equation  (17)  by  (g’  )‘  ,  we  find 


R’  =  R  +  R  (1  +  S  R  )’  , 
na  nt  np  mt  c 


(18) 


According  to  Schwantes’  results,  the  values  of  R  /R  were  always  of  the 
°  np  nt 

order  of  0.1  and  g  / g  was  alwavs  less  than  0,1,  so  that  it  is  negligible  for 
"x  mp 

most  practical  puposes. 

c)  Feedback  effect  with  correlation  between  i  and  i 
_ _ _ _ _ _  £  P 

Next  we  consider  correlation  between  i  and  i  ,  but  no  feedback.  Since 

c  p’ 

i  =  i  A  +  i  ,  the  spectrum  of  the  anode  current  is 
a  c  p’ 


i  -  =  i‘V+2ii  A  +  i ■ 
a  c  c  p  p 


(19) 


and  i  i  =  2qleq  cor,  f;  so  one  obtains 
c  p 


i  =  2q(I  V  +  2loq  cor.  \  +  I  )Af  .  (19a) 

a  eq  Cq 

c  Hp 

This  is  without  feedback.  If  we  consider  feedback  and  correlation,  we  have 

1,:  =  i,;''+2iT+r;  i*  =  - ^ -  (20) 

a  c  c  p  p  c  ,  +  SmtRc 
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and  we  can  get  the  final  formula 


=  2qAf 


r  i  *2 

eqc 

-(1+smtRc> 


+ 


2T  X 

eg  cor.  ,  T 

(1  +  g  R  )  eq 
&mt  c  ^p 


(20a) 


d)  Transformer  feedback  from  screen  grid  to  control  grid 

As  mentioned  before,  partition  1/f  noise  flows  from  screen  grid  to 
anode.  The  partition  1/f  noise  can  be  eliminated  by  connecting  a  transformer 
with  proper  turn  ratio  n  between  screen  grid  and  control  grid.  The  details 
are  given  in  C.J.  Hsieh's  Engineer's  thesis.  The  transformer  feedback  measure¬ 
ment  is  shown  in  Figure  31. 

e)  Experimental  setup  and  device  description 

The  circuit  of  the  measurement  system  is  shown  in  Figure  32.  The  noise 
devices  are  three  pentodes,  6AUGWC  (RCA),  6AG5  (RCA),  and  6CE6  (GE) .  Since  the 
pentode  is  a  high-gain  amplifier,  we  have  to  be  careful  with  the  input  terminal. 
Figure  32  shows  the  biasing  circuit  of  the  device  under  test.  Since  the  control 
grid  is  very  sensitive  to  60  Hz  induced  at  the  input,  a  low-pass  filter  is  used, 
so  as  to  reduce  the  60  Hz  noise  as  much  as  possible.  The  biasing  circuit  with 
low-pass  filter  consists  of  a  45V  DC  battery,  a  100 K  potentiometer,  and  a 
high-quality  paper  capacitor  (2  PF)  to  give  the  correct  bias  voltage  needed. 
Similarly,  we  choose  a  wire-wound  resistor  Rc  in  the  cathode  lead,  in  order  to 
reduce  the  termal  noise  as  much  as  possible.  Two  70V  batteries  give  the  cor¬ 
rect  screen  grid  voltage.  Six  70V  batteries  are  used  as  a  plate  power  supply. 
The  anode  current  I  is  monitored  by  an  HP  solid-state  digital  voltmeter. 

A  0.1  UF,  250  VW.V  coupling  capacitor  is  set  in  between  the  load  resistor  Rp 
and  the  Hewlett  Packard  3582A  spectrum  analyzer  to  prevent  the  spectrum 
analyzer  from  damaging  by  the  high-load  voltage.  The  HP  3582A  spectrum  analyzer 


featuring  a  dual-channel  Fast-Fourier  Transform  is  employed  essentially  as  a 
narrow-band,  tunable  filter. 
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Two  I-V  measurements  were  done  on  each  pentode.  The  DC  current- 
voltage  characteristics  of  the  pentode  6AU6WC  (RCA)  have  the  conventional 
form. 

The  whole  system  is  first  operated  without  cathode  feedback  resistor. 

The  variable  potentiometer  (see  Figure  32)  is  set  to  produce  different  biasing 
voltage  and  get  different  anode  current.  The  room  temperature  noise  current 
spectrum  (less  than  25  KHz)  for  three  anode  currents  is  shown  in  Figure  33 
[for  6AU6WC  (RCA)]. 

We  next  consider  the  measurement  with  cathode  feedback  by  inserting 
several  wire-wound  resistors  R  as  shown  in  Figure  32.  The  low-frequency  noise 
spectra  were  measured  for  several  resistors  (R^  =  300 ft,  500  ft,  700  ft ,  1  Kft, 

1.5  Kft,  2  Kft,  and  3  Kft)  up  to  3  Kft  as  shown  in  Figure  34  for  6AU6WC  (RCA).  The 
magnitude  decreases  with  increasing  resistance  up  to  about  3  Kft;  further  in¬ 
creasing  of  the  cathode  feedback  resistor  does  not  affect  the  magnitude  very 
much,  but  the  slope  becomes  less  than  1/f  beyond  100  Hz.  Figure  35  represents 

the  same  results  for  the  pentode  6CE6  (GE)  which  has  a  g  about  2.6  milimho 

and  feedback  resistor  increasing  from  500  ft  up  to  10  Kft.  Table  1  lists  the 

data  obtained  for  several  different  cathode  feedback  resistors.  The  noise 
spectral  density  is  obtained  from  the  result 


SI  7 


(1  +  g  R  ) 
mt  c 


(21) 


Experimental  results  are  compared  with  theoretical  values  for  the  pentode 

6AU6WC  (RCA)  at  100  Hz  and  1,000  Hz,  and  the  comparison  is  given  in  Table  1 

and  Table  2,  respectively.  The  comparison  of  the  curve  of  S^/S^  versus 

R  with  the  theoretical  plots  at  100  Hz  and  1,000  Hz  is  shown  in  Figures  3b 
c 

and  37,  respectively.  Ue  note  that  the  two  curves  almost  coincide.  Similarly, 
the  procedure  is  repeated  for  another  pentode  6CE6  (GE) .  Experimental  results 
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are  plotted  and  compared  with  the  theoretical  values  given  in  Table  3  and 
Table  4  and  are  shown  in  Figures  38  and  39.  As  mentioned  previously,  we  have 


sz  (f) 


sT  (f)  =  — 

a  (1 


- - Z  +  S  (f)  =  2qA f  f - ^ - r  +  b] 


(22) 


where  A  and  B  are  constants.  Equation  (22)  can  be  rewritten  in  terms  of 
equivalent  saturated  diode  current. 


Ieqa  = 


+  B  . 


(22a) 


(1+&mtRc>‘ 


Thus,  suppose  we  plot  theoretically  log  versus  log  (1  +  g  Rc)  >  the  curve 
is  identical  to  a  plot  of  log  (Ieqa)  versus  log  (1  +  g^R  ).  Therefore,  by 

T"“  o 

comparing  the  experimental  curve  of  log  (i^ )  versus  log  (1  +  g^R  )■*■  and 
log  (i^/ia^)  versus  log  (1+g^  R£)2  with  the  theoretical  curves  for  the  pen¬ 
tode  6AU6i vC  (RCA),  we  get  a  result  shown  in  Figures  40  and  41,  respectively. 
According  to  the  experimental  data,  we  can  calculate  the  A  and  B  values  as 
follows.  From  Equation  (22a)  and  assuming  we  have  (l  +  g^R^ “»1, 

so  that  B  -  Ieqp  and 

2cleqp  *  2  x  1.6  X  1(T19  X  teqp  =  3  x  IQ-21  , 


hence 


-19 

B  =  Ieqp  =  3  x  io  x  10  =  9.6  mA  . 


As  mentioned  before, 

X2  =  (I  /I  )2  =  (3.4/3.9)2  =  0.76  (1  *  3.4  mA,  I  *  3.9  mA) . 

3  C  3  C 

From  the  experimental  curve,  we  can  find  that  B  is  less  than  A.  So  it  is 
negligible  for  most  practical  calculation.  So 

A  =  IeqcV  =  Teqc  x  0.76 


2qleqc  x  o.76  =  S  =  1.224  x  10 
o 


-19 


and 


p 

L 
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I  =  (1.224  x  10-19)  /  (0.76  x  3.2  x  10-19)  =  503  mA 
eqc 

and 

A  =  I  =  503  mA  x  0.76  =  381  mA  . 

eqc 

If  Equations  (22)  and  (22a)  do  not  fit  too  well,  then  we  have  to  con¬ 
sider  the  possibility  of  correlation  between  i  and  i  .  Results  for  the  occur- 

c  p 

rence  of  this  are  found  in  Hsieh's  thesis.  See  also  Figure  42. 

Results  including  transformer  feedback  are  shown  in  Figure  43. 

f)  Effect  of  anode  voltage  on  partition  1/f  noise 

Now  we  observe  Figure  44;  it  has  feedback  but  no  transformer.  We 

measure  i?  at  the  higher  R  at  V  =  133.74  V  and  V  =  453.6  V.  In  each  case 
a  c  a  a 

we  do  ten  times  and  then  calculate  the  average  in  order  to  improve  the  accuracy. 

The  measurement  of  S^.  may  be  done  with  a  precision  of  1.2%  or  better. 

xa 

Here  we  measure  at  V  =  133.26  V,  V  =  133.74  V,  and  at  V*  =  133.43  V, 

8  a  g 

V'  =  453.6V.  This  hardly  changed  (from  5.78  mA  to  5.81  mA) ,  but  changed 

cl  C 

I  (from  3.97  mA  to  4.25  mA)  and  hence  1^  changed  from  1.81  mA  to  1.554  mA. 

a  g 

As  discussed  in  section  a)  for  the  partition  noise,  Handel  presents  the 
expression 

2\A 

S(f)  =  (l-r2)I-'  .  (23) 

f  8 

But  van  der  Ziel  suggested  that  the  partition  noise  nay  be  given  as 
2>r\/A,  A 

S(f)  =  - I  I  (23a) 

f  a  g 

so  that 

tS;(Wa 

[sp(f)]va 


(24) 


tlrv.uiv  that  a  factor  N  is  necessary  in  the  denominator  of  Kquation  (~) 


see 
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Equation  (25).  Van  der  Ziel  indicated  that  Equation  (7),  without  the  factor 

9 

is  off  by  a  factor  10  .  For  let 

SI  (f)  ’  T  V2  A?/"2  (27) 

P 

where  A  =  1.224  x  10  V  ,  A  9  =  1.224  x  10  ~*V„,  a  =  1/137  and  1  and  I. 
wd  a  z  a  z 

are  in  amperes.  Now  if  Vfl  =  V?  =  133.6  volts,  I=  =  5  mA,  I9  =  2  mA, 


Equation  (27)  yields 


Sj  (f)  =  (2.4  x  10  u/f)A?/H2  . 


Experimentally,  Hsieh  observed 

Sj.  (f)  =  (4  x  10_19/f)A2/H2  . 

P 

O 

This  yields  a  ratio  of  6  x  10  !  Thus  it  is  certain  that  a  factor  N  of 
9 

order  0.6  x  10  ,  as  in  Equation  (25),  must  occur.  This  N  must  be  the  number 

—  3  16 

of  electrons  in  the  system.  If  I  =  7  x  10  A,  then  I  /e  =  4.4  x  10 

c  c 

electrons/sec.  The  transit  time  is 

t,  =  2d/6  *  107VV  =0.48  x  10“8  sec.  (V  =  50V,  d  =  ]  cm), 

d  »  av  av 


N  =  (I  /e)T,  =  4.4  x  1016  x  o.48  x  icf8  =  2  *  108  . 
c  d 

g 

This  comes  very  close  to  the  experimental  value  involving  the  factor  6  x  iol  . 
Thus  (25)  is  the  correct  result  and  not  (27). 


Not_e:  It  should  be  noted  that  the  missing  of  large  factors  X,  as  indi¬ 

cated  here,  is  one  of  the  profound  problems  of  Handel's  theory.  The  same 
applies,  in  our  opinion,  to  Handel's  theory  for  emission  noise  involving  radio¬ 
active  decay  (see  Section  IV),  There  a  factor  of  N  *  500  would  be  required  in 
order  to  give  a  more  satisfactory  agreement  with  experiment.  The  apparent 
coherence  properties  of  scattering  particles  in  Handel's  theory  is  a  fact  that 
has  not  yet  adequately  been  solved.  This  remains  one  of  the  main  obstacles  of 
Handel's  theory. 
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C.  THEORETICAL  WORK 

I.  Remarks  on  quantum  1/f  noise  in  metal  films.  C.M.  Van  Vliet 
Hooge's  formula  reads 

SjCO/l2  =  a/fN  .  (1) 

In  nondegenerate  semiconductors  N  is  the  number  of  carriers  contributing  to 
the  noise  in  the  entire  device.  For  metal  films,  N  has  been  interpreted  as 
the  number  of  atoms  in  the  film,  multiplied  by  the  valency.  On  the  contrary, 
we  believe  that  N  is  to  be  the  number  of  electrons  near  the  Fermi  surface 
(a  slice  ~  2kT  around  the  Fermi  level) ,  since  the  other  carriers  cannot  con¬ 
tribute  to  mobility-fluctuation  noise.  An  intuitive  reasoning  was  presented 
in  Section  B-I.  A  more  complete  theory  based  on  the  Boltzmann  transport 
equation  is  being  developed.  A  preliminary  letter  to  the  Editor,  submitted 
to  Physica,  is  included  herewith. 
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Letter  to  the  Editor 


1/f  Noise  in  Mobility  Fluctuations  and 
the  Boltzmann  Equation 

A.  van  der  Ziel^  ,  C.M.  Van  Vliet^  , 

Dept,  of  Electrical  Eng. ,  U.  of  Florida,  Gainesville,  FL,  32611,  USA 

and 

R.J.J.  Zijlstra 

Natuurkundig  Laboratoriura  Ryksuniversiteit ,  Utrecht,  The  Netherlands 

and 

3) 

R.  Jindal  ' 

Dept,  of  Electrical  Eng.,  U.  of  Minnesota,  Minneapolis,  MN,  55455,  USA 
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Summary  J 

The  1/f  noise  in  mobility  and  diffusivity  are  discussed.  The  two  noise 
densities  are  related  due  to  an  instantaneous  Einstein  relation,  valid  for 
times  larger  than  the  collision  tine.  The  Boltzmann  equation  treatment  shows 
that  Hooge's  relation  for  the  collective  band  mobility  noise,  S^/(u)''  =a/fN, 
can  be  justified  for  a  nondegenerate  and  spatially  noncorrelated  electron  or 
ion  gas. 
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1.  Introduction 

Several  recent  papers  have  discussed  the  result  of  1/f  noise  for  the 
relative  mobility-fluctuation  spectrum  Sy(f)/(y)2  of  carriers  in  semi¬ 
conductors1  Some  papers  equate  it  to  a/f  and  others  to  a/fN;  here  a  is 
the  Hooge  parameter,  f  the  frequency,  and  N  the  number  of  carriers  in  the 
sample.  The  discrepancy  is  easily  resolved  by  distinguishing  between  the 

relative  1/f  mobility-fluctuation  spectrum  S  (f)/(y.)2  of  a  single  carrier 

^i  1 

and  the  1/f  mobility-fluctuation  spectrum  S^(f)/(y)2  for  an  ensemble  of 
carriers.  This  was  shown  recently  by  Hooge^.  The  pertinent  relations  are 
found  to  be 

S  ( f ) / (u . ) 2  =  a/f  (1.1) 

y .  i 

and 

S^(f)/(y)2  =  a/fN  .  (1.2) 

In  this  note  we  point  out  once  more  that  (1.1)  is  the  basic  relation 

and  that  (1.2)  follows  for  nondegenerate  semiconductors;  it  may  not  hold  for 

metals,  degenerate  semiconductors,  and  in  ionic  solutions  in  which  carriers 

7  8 

within  a  certain  distance  are  correlated  ’  .  Our  proof  for  the  validity  of 
(1.2)  in  nondegenerate  semiconductors  is  based  on  the  Boltzmann  equation. 

We  also  derive  an  Einstein  relation  for  the  fluctuations  and  find  the  corres¬ 
ponding  expression  for  the  relative  diffusivity  noise  SD(f)/(D)2. 


2 .  Boltzmann  equation  approach 

In  the  k-space  approach  we  start  from 


Jn(r,t)  =  -q 

where  J  (r,t)  is  the  electron  current  density  Z(k)  =  1/4  ;r2  is 
of  states,  and  v^  is  the  velocity  in  a  one-particle  state  |k); 
the  occupancy  of  a  state  k  in  the  vicinity  of  r,  averaged  over 


I 


d JkZ(k) f (r ,k, t)v 


(2.1) 

the  density 
f(r.k,t)  is 
many  collisions. 
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but  still  fluctuating  due  to  1/f  noise.  It  can  be  computed  from  the  mode 

occupancy  operator  average  < cj"  c  ,>  of  the  quantum  field,  compare 

k.  tc  t 

9  ~  ~ 

van  Vliet  et  al  . 

More  simply,  we  obtain  f  from  the  Boltzmann  equation,  assuming  that 

the  cross  section  for  scattering,  and  thus  the  relaxation  time  x  fluctuates, 

1  9 

either  due  to  the  fluctuating  phonon  population  ’  ,  or  due  to  quantum  1/f 
noise^.  Thus,  putting  f(r,k,t)  =  f^(r,k)  +  f3(r,k,t)  where  f^(r,k)  is  the 
local  equilibrium  distribution,  we  find  in  the  usual  way  in  the  relaxation 
time  approximation,  keeping  f°  in  the  streaming  terms  and  f1  in  the  collision 
integral, 


Jn(r,t)  =  q J  d3kZ (k) x(r,k, t) 

[(vkvk)  •  7f0(r,k)  -  (q2/h)vkE  •  Vkf°(r,k)]  . 


(2.2) 


Writing  still  ^f0  =  hvk3f°/3  <£^_,  this  yields  the  diffusion  part 

J  ■ . ,,  =  qV  •  (D  n)  [for  nondegenerate  semiconductors,  V  •  (D  n)  =  D  •  Vn, 
-n.diff  ~  -mi  h  ~  --n  =n  ~  ' 

since  D  does  not  depend  on  the  chemical  potential^]  and  the  drift  part 
=n 

J  ,  .  ,  =  qu  •  E  with 

~n, drift  ^n 


D  (r,t)  =  «x(r,k,  t)v.  v  » 

~  ~  -  K.~K 

\  ,  ,3  ioc  f  ® 

JJ  (r,t)  =  -  q  «x(r,k,t)  — r-p -  v  v.  » 

■  n  -  ~  3  c,-k  ~K~ic 

where  «  »  denotes  a  k-snace  average;  the  local, 
in  k-spacc  is  defined  by 


(2.3) 

(2.4) 


tone- 5cuC.tua.tuig  average 


«  1  ( r ,  t )  » 


n(r) 


/ 


'kZ(k) r(r,k,t)f °(r,k) 


(2.3) 


where  we  notice  the  occurrence  of  the  factor  n(r)  in  the  denominator. 


since  .  distribution  in  k-space  is  normalized  to  n: 


/ 


d3kZ(k)f° 


=  n(r)  , 


(2.6) 


(the  bar  denoting  a  time  average).  Restricting  ourselves  further  to 
scalar  D  and  p  and  nondegenerate  homogeneous  semiconductors,  for  which 
3  log  =  -1/kT,  eqs .  (2.3)  to  (2.5)  yield 


D 

n 


U 


n 


d3kZ(k)(v2/3)  T(r,k,t)f°(r,k)  , 

-  k  ~  ~  -  - 

d3kZ(k)(q/kT)(vk2/3)  r(r,k,t)f°(r,k)  . 


(2.7) 


(2.8) 


These  results  show  that: 

(a)  The  fluctuating  mobility,  as  found  by  the  Boltzmann  equation,  is 

of  the  form  (1/N)  Y.p  (k,t)  where  p  (k,t)  refers  to  the  mobility  of  the 
“K  n  •  n  ~ 

state  |k);  thus,  if  the  spectrum  for  un(k,t)  has  the  form  (1.1),  then  the 
spectrum  of  the  mobility  for  the  entire  conduction  band  has  the  form  (1.2) 
if  the  carrier  1/f  fluctuations  are  independent,  as  in  nondegenerate  semi¬ 
conductors.  For  clearly,  with  S  .  .  given  by  (1.1) 


_  _L  y  - 

p  ,  .  m2  i  w  (k) 

n,band  N"1  k  n  - 


N 


(2.9) 


A  more  accurate  computation  of  the  band  integrals  involved  shows  that  (2.9) 

1° 

must  be  modified  by  a  factor  of  3^/8,  cf,  Kleinpenning  and  Bell  “.  However, 
in  view  of  the  experimental  uncertainty  of  the  Hooge  parameter,  this 
correction  is  hardly  significant. 

(b)  There  is  an  instantaneous  Einstein  relation  for  the  fluctuating 

quantities  D  (t)  and  p  (t): 

n  n 


D  (t)  =  (kT/q)p  (t)  ; 
n  n 


(2.10) 
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Note  that  t  »  t  due  to  our  assumptions  of  a  vanishing  collision  integral 
for  f°(r,k,t)  =  f°(r,k).  We  thus  have  the  spectral  relation 

SD  (f)  =  (kT/q)2S  (f)  (2.11) 

n  Pn 


or 


sd  /(V2  ’  s„  /(v>2  •  (2-12) 

n  n 

13  14 

This  equation  was  conjectured  before  ’  but  never  proven. 

For  metals  and  for  ionic  solutions  above  certain  concentrations  a 
coherence  volume,  in  k-space  or  in  position  space^,  comes  into  play.  Then 


At(r,k,t)AT(r' ,k' ,t)  ^  F(r , k) 6 (r  -  r ' ) <$ (k  -  k ' ) .  Clearly  then,  though  (1.1) 
might  be  valid  for  the  individual  state  mobilities,  the  band  mobility  fluctua¬ 
tions  do  not  obey  (1.2).  A  -omputation  for  degenerate  semiconductors  and 
metals  will  be  published  elsewhere. 
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II.  1/f  noise,  due  to  Impurity  scattering.  Ganesh  Kousik 

A  calculation  is  made  of  the  mobility  (impurity  scattering  limited) 
fluctuations,  assuming  that  the  fluctuations  in  the  scattering  ion  section 
are  given  by  Handel's  theory  of  fundamental  1/f  noise. 

We  start  with 


co 

J_(t)  =  — f  e  Ww  ^  t(w,t)dw 

n  3m  V7  J 


and 


-7 - r  =  N  v 

x(t,w)  o 


/(I  -  cos  -1 


A)o^(t  ,w)dp. 


(1) 


(2) 


(within  the  relaxation  time  approximation) . 


=  Vw  /(1  -  003  •)A0o(t’w)d: 

— — —  =  N  2  v 2  [( 1 

Tlt  o  w  J 


cosv)(l  -  cos  :j')  <i0p(t,w)A0o,  (t,w)  >  dQdil'  (3) 

Assuming  that  scattering  in  different  directions  are  not  correlated  and 
finding  the  frequency  spectrum  on  both  sides  of  (3),  we  get 
S  (w) 


>  w  r 

- -  =  N  '  v  2  1(1- 

x4  O  W  J 


and 


So(^,Q',w)  =  <o(!!,w)>  ■iy^6(n-n*) 


cos  0)  (1  -  cos  O')  S  (n,a'  ,w)df?df2' 
o 


2.u\ 


(4) 


s  r 

—  =  N  2  v  2  [ 

t4  O  W  J 


(1  -  cos  0)  <o(P.,w)>2  ~  dn  . 


(5) 


Similarly,  assuming  that  scattering  of  electrons  with  different  energies  are 
not  correlated,  we  get 


CO 

S  =  ( — P~-\  f  e”2Ww3S  (w)dw 

U  \  3m"  V^)  { 


(6) 
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According  to  Handel’s  theory 


2aA  = 


16a 


(p  +mk  T/p  )  sin2  0/2 


(7) 


3™  c~ 


{mk^T/po  term  is  neglected},  9  is  the  scattering  angle,  and  p^  is  the  electron 
momentum.  Further, 


<a(fi,w)  >  = 


* 

2m  z 


-,4 


4iTett2  (4k^  sin2  6/2  +  k2) 


(8) 


for  a  screened  coulomb  potential,  k2  is  the  screening  parameter  which  for 
nondegenerate  semiconductors  is  the  Debye  screening  and  for  metals  is  given  by 
the  Thomas-Fermi  screening.  Impurity  scattering  is  assumed  elastic  (needed 
for  the  relaxation  time  approximation)  and  therefore  ]k-k'|  =  2 )  k  J  sin  9/2. 

Using  (5),  (6),  (7)  and  (8),  we  get 

8' 


-  / 


-2w 


_ _  J"  r  tan7 8 

■  -(>  -  ^  j]  u 


A  de 


dw 


(9) 


where  8'  =  tan  ^2k/K,  and  k2  =  2mt/-fi‘:  for  spherical  energy  surfaces. 


A’  = 


2-(16)4*ae4kB4T4 
27:im  '  c-  N  2  z4e6 


Equations  (1),  (2),  and  (8)  give 

c-o 

“«  '  B/ 


-W  D 

e  w 


o  £n  I  1  + 


dw 


(10) 


where 


;/2  .  ’/■■ 


B  = 


64  ^  v  T 

3  VnT  N  z  ‘  e  2 

f  n 
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We  now  have 


1  A' 
f 


/ 


~w 


w^dw 


U  o 


£n(l  +  4k2/ic2) 


A-ne- 

£kBT 


",2-t 


where  n  is  the  number  of  conduction  electrons  and  is  the  donor  concentra¬ 
tion.  e  the  dielectric  constant  is  11.9  c  . 

o 

At  temperatures  in  the  10°  -  70°K  range,  not  all  donors  are  ionized. 

To  find  n,  we  solve  the  neutrality  equation 

P  -  n  +  N  +  •-  0  (12) 

.  Therefore,  the  equation  for 

n  is 

n2  +  n^n  -  NQn^  0  (13) 

neglecting  terms  with  n.2  s;0.  Figure  45  shows  S  f/y2  vs.  T. 

t  u  n 

Mn 

In  metals  one  uses  the  Thomas-Fermi  screening 


where  N 


Vi 


D 


+  n 


and 


=  N  exp 
c 


m 


<2  =  4tto vi^Tt:-,)  (in  CCS  units) 
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where  the  density  of  states  at  the  Fermi  level. 


For  gold: 


22  cm 


_  .  „  L11I 

N  =  5.9  x  10 

£F  "  kBTF  =  kB  x  6*39  x  1C)4  er8s  • 

S 

yn 

and  spherical  energy  surface  (k^  =  2mc/h^)  have  been  used.  ■  —  vs.  T  for 

^n 

gold  are  shown  in  Figure  46. 

S  f 

Equation  (11)  implies  — yr  ,  the  noisiness,  is  not  a  strong  function  of 

P 

the  doping  density  in  semiconductors.  It  is  nearly  a  linear  function  of  T  in 
both  Si  and  gold,  in  fact  in  semiconductors  and  metals. 

The  noise  in  gold  (Fig.  45)  is  smaller,  probably  due  to  a  much  larger 
screening  constant  than  in  silicon.  Thus  the  electron  does  not  see  the  ion 
potential  as  strongly  as  in  Si. 

For  a  purely  coulomb  potential,  we  find  x  -►  0  and  therefore 


S  -  0  ; 
Pn 


also, 

s 

J  r'*  ...  dw 

P  -*•<>,  but  7  "*■  ” 

I 

,  [ 

I  T  .  . .  dw  ‘ 

III.  Noise  Out  of  Equilibrium 
by 

Carolyn  M.  Van  Vliet 


Department  of  Electrical  Engineering 
University  of  Florida,  Gainesville,  32611 

and 

Centre  de  Recherches  Mathematiques 
Universite  de  Montreal,  Montreal,  Quebec  H3C3J7,  Canada 


Summarv 
— . , 

We  discuss  six  properties  of  stochastic  processes  which  are  generally 
valid  only  in  thermal  equilibrium  systems.  In  Part  A  of  this  paper  we  discuss 
the  first  and  most  fundamental  property,  viz,  the  fluctuation-dissipation 
theorem.  Kubo's  and  Van  Vliet’s  quantum  statistical  mechanical  derivations 
are  reviewed.  Next  we  discuss  Gupta's  thermodynamic  theorem,  which  may  hold 
for  a  class  of  nonequilibrium  systems.  In  Part  B  we  discuss  the  other  proper¬ 
ties  which  mainly  center  around  the  principles  of  mesoscopic  reversibility  and 
detailed  balance,  and  the  syrametrization  property  of  the  two  terms  in  the 
generalized  second  moment  relaxation  theorem  (also  called  generalized  g-r 
theorem,  generalized  Einstein  relation,  or  A-theorem).  We  then  give  two  ex¬ 
amples  for  which  this  syrametrization  property  does  not  hold.  First  we  review 
Van  Vliet ’s  model  of  nonequilibrium  optical  pumping  in  photoconductors  with 
traps.  Finally,  we  review  Tremblay's  et  al.  derivation  of  Brillouin  light 
scattering  as  an  example  of  a  nonequil ibrium  hydrodynamic  steady  state. 


1 .  Thermal  equi Librium  properties 

There  are  a  number  of  properties  which  render  the  description  of  noise 
in  thermal  equilibrium  considerably  simpler  than  in  a  nonequilibrium,  or  driven 
equilibrium,  steady  state.  The  main  features  for  fluctuation  processes  are  as 
follows. 

"Invited  paper  presented  at  the  7th  International  Symposium  on  Noise  in 
Physical  Systems,  Montpellier,  France,  May,  1983. 
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1)  Validity  of  the  fluctuation-dissipation  theorem 

2)  The  principle  of  mesoscopic  reversibility 

3)  The  principle  of  detailed  balance 

4)  The  symmetrization  property  of  the  "generalized  second  moment 
relaxation  theorem" 

5)  The  equilibrium  Einstein  relation  between  generalized  mobility  and 
dif fusivity 

6)  The  Onsager  relation  for  reciprocal  flow  problems. 

The  derivation  of  (1)  will  be  briefly  given  below.  The  fluctuation- 
dissipation  theorem  is  a  genuine  consequence  of  the  existence  of  the  micro- 
canonical  or  canonical  ensemble  to  describe  equilibrium  properties.  It  can  be 
derived  from  a  purely  thermodynamic  point  of  view  (Nyquist  1928,  Callen  and 
Greene  1952)  or  from  a  quantum  statistical  mechanical  approach  (Kubo  1957). 

The  other  five  properties  are  somewhat  more  subtle  and  have  a  less 
general  form  in  that  we  must  distinguish  between  a-variables ,  which  are  in¬ 
trinsically  time  reversible,  such  as  the  position,  and  ^-variables,  which  are 
odd,  i.e.,  change  sign  under  time  reversal,  such  as  the  velocity  (Casimir  1945). 
The  details  of  this  necessary  distinction  are  often  overlooked,  and  a  number  of 
errors  occur  even  in  recent  literature.  A  much  more  careful  consideration  than 
hitherto  given  is  necessary,  and  errors  occur  both  in  Lax  (1960)  and  the  author's 
papers  (1965)  in  this  respect.  Since  the  space  of  this  conference  paper  is 
limited,  we  shall  make  no  attempt  to  rectify  these  errors  hereby,  restricting 
ourselves  to  a-variables.  Some  remarks  on  the  other  cases  will  be  made  in  some 
instances.  The  five  properties  2)  to  6)  all  hang  together  in  a  unique  way.  One 
can  show  that  the  validity  of  any  one  of  the  five  implies  the  validity  of  the 
others.  Again,  the  writing  of  a  comprehensive  survey  paper  on  this  topic  would 


be  extremely  useful,  but  space  is  lacking  here  since  we  have  been  asked  to  write 
on  nonequilibrium. 
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This  paper  will  therefore  be  divided  as  follows.  In  Section  2  we 
review  the  standard  derivation  due  to  Kubo,  and  its  recent  modification  by 
Van  Vliet  (1978,  1979)  of  the  fluctuation  dissipation  theory.  We  also  discuss 
how  a  true  general  extension  to  nonequilibrium  systems  should  be  obtained. 
Lacking  such  a  general  extension  at  this  moment,  we  review  in  Section  3 
Gupta's  thermodynamic  attempt  to  obtain  a  meaningful  extension  for  electric 
nonlinear  circuits  under  a  set  of  restricted  but  useful  conditions  (Gupta 
1978,  1982). 

In  the  rest  of  the  paper  we  discuss  the  extensions  for  nonequilibrium 
when  any  of  the  properties  2)  to  6)  does  not  hold.  This  is  relatively  easier, 
since  full  solutions  for  Markovian  nonequilibrium  systems  have  existed  since  the 
initial  work  on  optically  excited  generation-recombination  statistics  in 
semiconductors  by  Van  Vliet  and  Blok  in  1956  (a  & b) .  The  first  true  non¬ 
equilibrium  example  involving  the  "generalized  nonequilibrium  second  order 
moment  relaxation  theorem"  (also  called  generalized  Einstein  relation  by  Lax 
1960,  generalized  g-r  theorem  by  Van  Vliet  and  Fassett  1965,  and  ft- theorem  by 
Van  Vliet  1971)  was  given  by  Van  Vliet  in  1964.  This  is  discussed  in  Sec¬ 
tions  4  and  5. 

In  recent  years  renewed  interest  in  the  nonequilibrium  problem  arose 
from  light  scattering  in  the  presence  of  a  small  temperature  gradient,  causing 
nonsymmetric  Brillouin  peaks.  This  problem  was  treated  by  a  number  of  authors 
(Procaccia,  Ronis  and  Oppenheim  1979,  Kirkpatrick,  Cohen  and  Dorfmann  1979, 
Tremblay,  Siggia  and  Arai  1980).  We  will  here  present  the  treatment  of  the 
latter  paper  (Section  6).  It  will  be  shown  that  their  example  is  mother 
straightforward  nonequilibrium  application  of  the  "generalized  nonequilibrium 
second  order  moment  re  1 axation  theorenT  as  originally  set  forth  by  the  author 
and  co’workers  (1954,  1965,  1971). 
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Part  A.  The  fluctuation-dissipation  theorem  and  extensions 
2.  Quantum  statistical  derivations 
Kubo  considers  systems  with  a  Hamiltonian 


H 


total 


=  H  -  AF( t) 


(2.1) 


Here  H  is  the  Hamiltonian  of  the  system  proper  and  -AF(t)  the  coupling  to  an 

external  field.  For  electrical  conductors,  F(t)  -*■  qE(t)  and  A  =  ]T  r.*  where 

i 

r^  are  the  position  operators  of  all  of  the  particles.  The  Von  Neumann 
equation  for  the  density  operator  now  reads 


-§J  +  (£|  [H,p]  =  (i)  u(t)F(t)  [A,p] 


(2.2) 


where  u(t)  indicates  that  the  field  was  switched  on  at  t  =  0.  If  we  are  now 
near  a  thermal  equilibrium  state,  we  can  substitute  p -+ p  in  the  r.h.s.  of 

— on  — pH 

(2.2)  where  p  is  the  canonical  density  operator  p  =  e  /Tr  e  ,  S  =  1/kT. 
eq  eq 

The  equation  (2.2)  can  easily  be  solved  for  the  average  response  of  a  flow 
quantity  B  where  B  is  an  operator  in  the  system.  One  finds  in  a  straight¬ 
forward  way  (Van  Vliet  1978) 

t 

<  2>B(t)  >  =  Tr  p  ( t ) B  -  Tr  ,  ^B  =  J  dT  •  t  (t  -  x)  F(t)  .  (2.3) 


di  .*  *A(t  -  t)F(t)  . 


Here  is  the  response  function  given  bv  the  commutator 

dA 

HaM  .-^TrMA,S(t))r,oql  (2.0 

where  the  time  dependence  of  p (Sclirodinger  picture)  was  transferred  to  the 
time  dependence  of  B  (Heisenberg  picture): 


•  .  i/t  S  iHt/-h  S  -i'rit/tf 

B(t)  =  e  B  =  e  Be 


(2.5) 


p  »« 

where  Bl>  is  the  Schrodinger  operator  (the  superscript  S  will  henceforth  be 

tted);/  is  the  Liouville  superoperator  acting  in  the  Licxiville  space  that 


omi 


contains  all  operators  A,  B,  ... 


In  particular,  let  §  be  the  current 
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J  =  B  =  r.  =  ^  A,  where  fi  is  the  sample  volume.  If  we  take  the  Founer- 
~  ~  SI  *!<  ~1  fi  ~ 

1 

Laplace  transform  of  (2.3),  then  we  have,  denoting  by  -  the  transform. 


j(ico)  =  q£(q/fl)A,A(ilj),g(lw) 


so  that  the  conductivity  becomes 


(2.6) 


Qo 

^(iuj)  =atr£if  dt 


e"Lujt  Tr  {  [A,A(t)]p 


eq 


CO 

/ 


-iwt 

=  SI  J  dt  6  ■  >  -  1  Tr  {  [J,J(t)]p  }  . 

'flu)  ~  ~  eq 


(2.7) 


We  have  here  the  beginning  of  Nyquist’s  theorem:  the  conductivity  is  related 
to  the  correlation  expression, 


<[J,J(t)]>  =  <  JJ(t)>  -  <  J(t)J> 


(2.8) 


where  J  J  =  J(0). 

Now  to  obtain  the  correlation  function  proper,  it  must  be  defined  as 
an  ant ieommutator ,  since  the  product  JJ(t)  must  be  symmetrized: 


( t )  =  <  [J,J(t)J  +  >  . 


(2.9) 


It  is  now  possible  to  relate  the  commutator  to  the  anticommutator  with  some 
trans format  ions  involving  contour  integration  for  the  Fourier  transforms.  It 
is  then  possible  to  express  the  full  (two-sided)  Fourier  transform  of  the 
correlation  anticommutator  into  the  single-sided  Fourier  transform  given  by 
(2.7).  The  result  is  the  spectral  expression 

S.  (...)  =  4£(w,T)  {(o’  (<J]s  +  i[a"  (u))]3}  (2.10) 

J  J.,  vp  vp 

where  v  and  p  denote  either  two  different  (but  correlated)  currents  or 
Cartesian  current  components,  o  =  o’  +  in",  and 


s  1 


a  1 


vp 


—  [o  +o  ],a  =  -r  [ a 

2  v«  pv  vp  2 


-  o  J 
vp  pv 
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further,  £(u,T)  =  (-W2)  coth  (31W2)  is  the  mean  energy  of  a  harmonic 
oscillator  of  frequency  mode  w.  Eq.  (2.10)  is  the  fluctuation-dissipation 
theorem  in  all  its  glory. 

We  have  criticized  Kubo's  derivation  in  our  1978  paper,  since  nowhere 
was  a  form  for  the  Hamiltonian  H,  commensurate  with  dissipation,  introduced! 

In  particular,  the  Heisenberg  operators  do  not  represent  the  required  approach 
to  equilibrium.  As  we  indicated  in  our  paper,  a  partitioning  of  the  Hamiltonian 
is  essential;  i.e.,  instead  of  (2.1)  we  write 

“total  =  H°  +  XV  “  AF(t)  •  (2.11) 

Here  H°  is  the  largest  Hamiltonian  that  can  be  diagonalized  for  the  many-body 
system,  e.g.,  an  electron-phonon  system;  XV  represents  the  interactions  which 
randomize  the  energy  over  the  states  of  H°.  In  the  Van  Hove  limit  (A-*0,t 
large.  At  finite)  and  large  system  limit,  such  randomization  leads  to  irrever¬ 
sibility  (The  Poincare  cycle  becomes  "off  limits"  for  observation  in  such  a 
system).  In  our  second  paper  on  linear  response  theory  we  indicated  that  (2.2) 

now  is  carried  over  in  an  irreversible  master  equation,  which  for  the  reduced 
£ 

p  (i.e.,  after  the  Van  Hove  limit)  now  reads 
,  R  R 

+  (A.  +  i  <£°)e  ( t)  =  F(t)  I  times  a  functional  of  p(t)  and  A].  (2.12) 

vC  Cl 

Here  A^  is  the  master  operator 

AdK  =  ~  £  j  V  >  <  V  !  W  ,  l  <  V '  i  K |  y '  >  —  <  y  |  K I  Y  >  ] ;  (2.13) 

ay' 

where  {|y>}  is  the  set  of  nany-bodv  states  of  H°,  is  the  interaction 

Liouville  operator  K  =  4)  ^[H°,K],  and  W  , ,  is  the  transition  probability 

i  V 

according  to  the  golden  rule,  W_,  =  (2  /«fi)  !  <  y  j  V  |  >’  >  |  '&  ( , )  . 

Equation  (2„12)  can  be  used  similarly  to  (2.2)  to  obtain  the  transport  co¬ 
efficients  a,  and  others,  to  be  expressed  in  correlation  functions  involving 
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R  R 

the  reduced  operators  J  (t)J(O),  where  now  J  (t)  =  exp[  -  (Aj  +  i  £?)t] J(0) . 

d 

The  main  problem  with  (2.12)  is  that  the  r.h.s.  has  so  far  only  been  evaluated 
for  a  near-thermal  equilibrium  state.  For  a  driven  equilibrium,  the  Hamil¬ 
tonian  (2.11)  should  be  further  partitioned: 

Htotal  =  H°  +  xv  “  AF  ~  Af(t)  ~  Fa(t)  •  (2.14) 

Here  AF  would  be  the  steady  driven  equilibrium  coupling  to  the  sustaining 
field  F  ,  and  f(t)  would  be  a  small  time-dependent  field  perturbation.  Could 
Equation  (2,12)  be  reformulated  for  this  Hamiltonian  with  only  the  part 
-Af(t)  occurring  in  the  r.h.s.,  then  we  could  obtain  a  steady-state  fluctuation- 
dissipation  theorem.  Such  a  program  has  not  yet  been  carried  out.  One  state¬ 
ment  can  be  easily  deduced  from  such  a  theory,  however.  We  note  that  in  the 
above  the  two  Hamiltonians  IV  and  AF  are  on  an  equal  footing.  The  noise  due 
to  both  is  a  function  of  the  small  signal  transport  coefficient  caused  by  f(t). 
Since  traditionally  W  causes  thermal  dissipative  noise  and  AF  causes  current 
or  shot  noise,  it  is  clear,  then,  that  in  a  generalized  nonequilibrium  fluctuation 
dissipation  theorem  we  wiLl  no  longer  be  able  to  separate  "thermal"  and  "shot 
noise";  both  will  be  contained  in  the  new  fluctuation-dissipation  theorem.  This 
is  born  out  by  the  theory  due  to  Gupta,  reviewed  below. 

3 .  Gupta's  thermodynamic  result 

Since  F  is  the  driving  parameter  which  is  fixed,  typically  a  large  d.c. 
field,  the  entropy  will  contain  intensive  variables,  as  foreseen  by  Gibbs 
(1902),  a  fact  overlooked  bv  Callen  in  his  well-known  monograph  (1962).  Thus, 
properly  the  entropy  is  a  Massieu-typo  Legendre  transform  of  Callen's  entropy. 

From  (2.14)  we  arrive  at  a  Gibbs  entropy  production 

Td  S  -  d  £  -  Ad  F  (3.1) 

4 

where  TdS  represents  the  dissipation  caused  by  XV.  This  form  satisfies 
Gupta's  thermodynamic  theory.  His  basic  assumptions  are  noted  in  his  1Q82 
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paper.  The  most  important  are: 

a)  The  entropy  depends,  besides  on  £,  on  another  variable  F,  which 
can  pass  through  the  port  of  the  system  and  which  is  fixed. 

b)  The  system  is  purely  resistive,  so  no  energy  is  stored  in  AdF.  We 

t 

should  therefore  be  able  to  write  A  as  a  function  of  the  instantaneous  values 
of  F:  A  =  CL<  F>  +  C0  <  F2  >  ,  etc. 

c)  The  system  has  no  quantum  correction.  This  is,  in  our  case,  born 
out  by  the  use  of  the  A  operator. 

For  the  present  case  of  electrical  conduction,  A  can  again  be  inter¬ 
preted  as  (f2/q)J  =  IL/q  where  I  is  the  current  and  F  =  qE  =  Vq/L.  Thus 
AdF  =  IV.  Gupta's  theorem  for  the  noise  now  reads 

S__(u>)  =  4kT  P  /  ,  (3.2) 

FF  ex 

or 

sWCj)  =  4kT  Pex  7  1 '  »  (3‘3) 

where  P  ^  is  the  excess  power  dissipation  caused  by  the  presence  of  a 
periodic  zero  average  small  signal  a  (or  i)  superimposed  on  the  driven 
steady  state.  For  a  nonlinear  resistor  with  V  =  V(I) 


Define  now 


Then  (3.6)  reads 

Sn<w)  =  4kT  [g(IT)  -  0IT] 


(3.7) 


(3.8) 


where  we  replaced  by  the  more  useful  symbol  1^  (for  total  average 
current) . 

Application:  Schottky  barrier  diodes  or  p-n  junctions  at  low  fre¬ 
quencies  (no  energy  storage  in  extra  term).  Then 


I  =  Il(e«V/kT  -  1) 


(3.9) 


g(IT)  =  (dI/dV)T  =  (q/kT)(IT  +  IL)  (3.10) 

SIT  «  (q/2kT)Ix. 

Thus  (3.8)  reads 

Su(cj)  =  2q  [  (IT  +  Ix)  +  IXJ  ,  (3.11) 

indicating  shot  noise  of  f orward  and  reverse  current.  Clearly,  this  com¬ 
prises  the  thermal  noise,  for  at  zero  bias 

SI1(~) I V=0  ”  4qIl  =  4kTglv=0  *  (3.12) 

As  we  noted  before,  no  distinction  can  be  made  between  thermal  noise  "proper" 
and  shot  noise. 

Whereas  Gupta's  theory  is  highly  specialized,  it  has  opened  the  way 


to  show  that  steady  state  dissipation  fluctuation  theorems  do  exist,  contrary 
to  previous  pessimistic  claims  (for  a  survey  see  van  Kampen  1963).  Gupta’s 
result  (3.6)  or  (3.8)  is  not  a  trivial  result.  Others  have  tried  to  generalize 
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the  fluctuation-dissipation-theorem  on  more  heuristic  grounds.  Van  der  Ziel 
(1973)  considered  the  following  possibilities  for  Nyquist's  theorem  in  non¬ 
linear  devices: 

a)  <vn:(f)>  =  4kT  B(V/I):  <i^(f)>  =  4kT  B(V/l)  |y  (f)  | 2 

b)  <in2(f^)>  s  4kT  B(I/V)j  <vn2  (f)>  =  4kT  B(I/V)  /  |y(f)  | 2 

c)  <vn°(f)>  =  4kT  B(dV/dI);  <i^(f)>  -  4kT  B  (dV/dl)  |  y  (f )  |  2 

d)  <  (f)>  =  4kT  B(dI/dV);  <  (f)>  =  4kT  B(dI/dV)  /  (y  ( f )  | ' 

e)  none  of  the  above 

In  view  of  Gupta's  result,  possibility  e)  is  the  answer  to  this  multiple  choice 
question,  even  though  it  can  be  shown  that  some  devices  satisfy  some  of  the 
other  possibilities  (which  then  must  be  equivalent  to  (3.6)  or  (3.8)). 


Part  B.  The  generalized  nonequilibrium  second  moment  relaxation  theorem 

(generalized  g-r  theorem,  generalized  Einstein  relation,  or  A-theorem) 


4.  The  mesoscopic  Markov  process 

In  statistical  mechanics  the  variables  are  the  microscopic  variables 
p^  and  q^  or  their  quantum  operators.  In  stochastic  processes  we  deal,  however, 
with  much  more  coarse-grained  variables,  which  arc  averages  over  large  ranges 
of  quantum  states,  even  though  the  fluctuations  in  these  variables  cannot  be 
seen  until  after  suitable  amplification.  Brownian  motion  and  fluctuations  in 
carrier  populations  (conduction  band,  traps)  are  typical  examples.  We  call 
these  variables  mesoscopic  (term  of  van  Kampen  1962).  They  can  be  pictured  in 
"a-space"  (De  Greet  and  Mazur,  1962). 

From  the  microscopic  master  equation,  see  (2.12),  one  can  derive  a 
mesoscopic  master  equation.  The  equation  roads 


.  J  dViP(.v,t|i,’)Q.1.1.1  -  p(.-.,th 


(4.1) 
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where  Q  „  is  the  transition  probability  per  unit  time  for  a  change  a"-* a. 
Assuming  that  a  is  a  set  of  variables  on  (0,°°),  one  can  Laplace  transform  the 
master  equation,  and  one  easily  finds  for  the  characteristic  function  (Van  Vliet 
1983) 


-s  *a(t) 


at  <e 


-s  *a(t) 


•tif  a”:  Fn(a(t))>a, 


(4.2) 


where  F  is  the  nth  order  Fokker-Planck  tensor 
n 

F  (a")  =  f  (a  -  a") nQ  „  dSa  =  <  a(At)  -  a"  >n„  /  At  . 
n  j  a^a  ~  ~  a 


(4.3) 


Differentiating  now  repeatedly  with  respect  to  s,  and  setting  s=0,  one  finds 
the  moment  equations.  The  first  moment  equation  is  the  phenomenological  equa- 


<Aa(t)  >  ,  =  -M<  -a(t)  >  , 

(sub  a '  means  a  conditional  average  in  an  ensemble  with  a(0)  =  a'  fixed); 
M  is  the  phenomenological  relaxation  matrix;  it  is  defined  bv 


(4.4) 


[ij  "  "  3aV  £  (ai"  ai)Qa"ala"  =  a0  * 
J  a  ~  ~J  ~ 


(4.5) 


where  a^  e  <a>,  the  unconditional  stationary  average.  By  differentiating 
twice  to  s,  setting  s  =  0  and  letting  t-*-00,  one  finds  from  (4.2)  the  generalized 
nonequilibrium  second  moment  relaxation  theorem; 


<  lala  >  N  +  M <  AaAa  >  =  B(a°) 


(4.6) 


where  B  is  the  second  order  Fokker  Planck  moment;  M  is  the  transpose,  la  is 
a  column  matrix  and  A a  a  row  matrix. 

The  theorem  was  first  derived  for  generation-recombination  noise  in  multi¬ 
level  semiconductors  or  photoconductors  (note  that  the  theorem  holds  out  of 
equilibrium  as  in  a  photoconductor!);  in  this  case  a  represents  the  various 
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populations  of  carriers,  a  =  (n^n2  •••  ns}c  The  theorem  was  called  the 
generalized  g-r  theorem  (Van  Vliet-Blok  1956a) .  For  three-dimensional  Brownian 
motion  the  theorem  takes  a  special  form  (M.C.  Wang  and  Uhlenbeck  1943).  The 
phenomenological  equation  for  this  type  is  dv/dt  =  -gv.  Thus  M+B.  The  second- 
order  Fokker-Planck  moment  is  easily  found  to  be  B_  =  2Dg2I,  where  I  is  the 
unit  tensor,  D  is  the  diffusivity.  Further,  <AvAv>  =  (kT/m)I.  Thus  (4.6) 
leads  to  mD  =  kT/g,  which  is  the  Einstein  relation  for  Brownian  motion.  (For 
carrier  motion  in  a  simple  semiconductor  model  g  =  1/t,  y  =  qt/m,  so  the  above 
leads  to  qD  =  pkT.)  For  this  reason  some  authors  (Lax  1960)  have  called  (4.6) 
the  generalized  Einstein  relation.  When  dealing  with  transport  processes, 

<AaAa>  is  to  be  replaced  by  the  covariance  kernel  <  Aa(r) Aa(r ' )  >  and  the 
corresponding  theorem  has  been  called  the  A-theorem  (Van  Vliet  1971).  To  us, 
at  present,  the  name  "generalized  nonequilibrium  second  moment  relaxation 
theorem"  seems  most  appropriate. 

We  now  state  the  equilibrium  properties  2)  -  6)  mentioned  in  the  intro¬ 
duction. 

a)  Microscopic  reversibility  W  ,  =  W  ,  follows  from  the  golden  rule 

YY  Y  Y 

in  quantum  mechanics.  The  mesoscopic  transition  probabilities  will  be  denoted 
by  Q  ,.  Let  now  x(a)  da  be  the  number  of  quantum  states  |y>  when  a  varies 

between  a  and  a  +  da.  We  then  clearly  have,  since  £  -*•  I  \(a)  da, 

Y  J 

x(a")Q  ..  =  x(a)Q  „  .  (4.7) 

cl  d  del 

This  is  the  principle  of  mesoscopic  reversibility.  In  the  microcanonical  or 
canonical  equilibrium  ensemble  the  probability  W(a)da  is  proportional  to  the 
number  of  accessible  quantum  states  in  da,  ice.,  x(a)da.  Thus  (4.7)  leads  to 
the  property  of  "mesoscopic  reversibility" 

W (a")Q  „  =  W(a)c)  „  .  (4.8) 

€  X  ci  t  1 
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b)  The  above  leads  to  detailed  balance  (Van  Vliet  1964).  Consider 
that  a  labels  the  occupancy  of  a  set  of  quantum  levels,  i.e.,  let  a  =  m,  n, 
or  k,  which  represent  discrete  population  vectors.  In  particular,  we  consider 
transitions  k  ++■  m  such  that 

{k,  ...k.,  k.  ...k  }  -*-►  fra.  . . .  m.  ,  m .  ...m  }  =  {k,  . . .  k.  +  1,  k .  -  1,  . . .  k  } , 
X  i  j  s  i  j  s  ^  i  i  j  s' 

(4.9) 

Let  p..  be  the  transition  rate  from  levels  .  to  levels  (governed  by  mass 

ij  t  j 

action  laws  or  similar  rates).  Then 


^km  =  Pji(^  5  Qmk  =  Pij(^  ' 


(4.10) 


Eq.  (4.8)  reads  in  this  notation 


-  "tv  • 


(4.11) 


We  sum  this  result  over  all  m.  Then  from  (4.9)  and  (4,10) 


Z  •••  Z  Z  **■  Z  w(-)pji(~} 


k  k.  k.  ,  k  n 

1=0  i=0  j=l  s=0 


I  -  I  I  -I  U (m) p^j  (m) 


ni.  a  m .  t  m .  in 

1=0  i  =  l  j=Q  s-0 


(4.12) 


or  also 


<  Pjt(k)  >  =  <  Pii  (j)  >  . 


(4.13) 


The  variables  k  and  m  are  dummy  variables;  thus  (4.13)  simply  says 


<P..>  =  <p..>.  This  is  detailed  balance  for  the  rates  between  anv  two 
Ji  ij 

sets  of  quantum  states  in  the  system. 

c)  From  Bayes'  theorem  one  easily  shows  that  (4.8)  leads  to  the 


symmetry  of  the  pair  correlation  function  for  small  intervals  At 


W  (a, At;  a",0)  =  W.,(a ",At;  a,0) 


(4.14) 
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With  the  linearized  phenomenological  equation  this  leads  to 

<a(At)a(0)>  =  <aa>  -  M<AaAa>At  (4.15) 

and 

<  a(0)a(At)  >  =  <aa>  -  <  AaAa  >'!At  .  (4.16) 

When  in  (4.14)  we  multiply  by  aa"  and  sum  over  both  variables,  we  easily  see 
that  the  S-.h.  sides  of  Eqs.  (4.15)  and  (4.16)  are  equal.  Thus ,  in  thermal 
equilibrium  the  r.h.  sides  of  (4.15)  and  (4.16)  must  also  be  equal,  which  is 
the  symmetrization  postulate  of  the  "generalized  second  moment  relaxation 
theorem";  in  (4.6)  both  sides  are  now  equal 

<AaAa>M  =  M<AaAa>  .  (4.17) 

This  property  gives  a  tremendous  simplification  for  thermal  equilibrium 
processes;  for  now  (4.6)  can  at  once  be  solved 

M  <  AaAa  >  =  -y  B(a°)  (4.18) 

or 

<  An  la  >  =  4"  M"1  B(a°)  .  (4.19) 

Since  (4.13)  is  for  Brownian  motion,  nothing  but  the  Einstein  relation  between 
mobility  and  diffusivitv,  it  is  clear  that  some  authors  (Lax)  refer  to  the 
full  theorem  (4.6)  as  the  nonequilibrium  generalized  Einstein  relation.  Well, 
whatever  the  name  of  Eq.  (4.6),  we  do  note  that  (4.6)  is  a  genuine  nonequilibrium 
result,  while  (4.18)  or  (4.19)  is  an  equilibrium  result.  Thus,  the  properties 

<  Aa.'.a  >  M  =  M<  AaAa  >  (4.20a) 

<  AaAa >  M  *  >1<  AaAa >  (4.20b) 

delineate  clearly  between  equilibrium  and  nonequilibrium  behavior. 

We  still  mention  the  analog  for  transport  processes  (Van  Vliet  1971). 
Suppose  we  have  a  stochastic  variable  depending  on  a  (vector)  parameter  y, 
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which  is  continuous  in  DC-'",00)  .  The  stochastic  process  a(y,t)  is  then  in¬ 
finite  dimensional,  though  we  assume  it  to  be  still  Markovian.  The  Langevin 
equation  is  of  the  matrix  operator  form 
3a(y,t) 

— — - +  Aya(y,t)  =  5(y,t)  .  (4.21) 

where  A  is  some  (matrix)  integral  or  differential  operator  describing  the 
transport  process.  Let  the  covariance  matrix  kernel  be 

JLCy.y')  =  <  Aa(y)Aa(v ')  >  .  (4.22) 

The  analog  of  the  matrix  theorem  (4.6)  is  now  the  A- theorem 

Ayr(y,y’)  +r(y,y’)Ay^  = s,- (y ,y’)  (4.23) 

t  r 

where  A  acts  from  the  right  on  T.  This  theorem  was  proven  using  Hilbert- 
space  methods  by  the  author.  Notice  that  S_  is  the  white  spectral  density  of 
the  Langevin  sources.  Again,  when  both  terms  on  the  f  .h.s.  of  the  .'.-theorem 
(4.23)  are  unequal,  we  deal  with  a  true  nonequilibrium  state. 

In  the  next  two  sections  we  give  both  a  matrix  and  a  transport  example 
of  a  true  nonequilibrium  situation,  which,  moreover,  has  been  experimentally 
verified. 


5 .  Optical  pumping  in  a  nonequilibrium  steady  state 

In  1964  we  considered  a  model  of  a  phc  t oconductor  in  which  free  electron 
hole  pairs  were  created  by  optical  absorption,  whereas  all  of  the  recombination 
was  occurring  via  intermediate  states,  see  Figure  1. 


Figure  1 


~7  //  '//>■.  \y'~; 


6  and  k  are  capture  constants.  Clearly,  there  is  no  detailed  balance,  so  we 


have  a  genuine  nonequilibrium  state.  The  M  matrix  follows  from  the  linearized 
mass  action  type  kinetic  equations  for  dAn/dt  and  dAi/dt  (Ap  =  Ai  +  An  is  a 
dependent  variable).  One  obtains 


One  can  now  solve  the  nonequilibrium  moment  relation  (4.6),  expressing  i  into 

n  via  (5.3).  Since  (4.6)  is  homogeneous  in  6/k,  the  relative  variance 
o 

<An:>/n^  can  be  computed  as  a  function  of  az(=«L)  for  various  values  of  the 
ratio  y  =  k/6.  The  results  are  given  in  Figure  2  below. 
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The  pecularity  of  this  model  is  that  <An2>/nQ  may  attain  values  >>1,  as 
are  indeed  experimentally  observed  in  cadmium  sulfide.  For  a  thermal  equili¬ 
brium  situation,  on  the  contrary,  we  always  have  < An2  >/nQ  <1,  according  to 
ermi-Dirac  statistics.  Thus,  the  model  of  this  section  represents  a  true  non¬ 
equilibrium  situation,  in  which  familiar  concepts  like  Fermi-Dirac  statistics 
are  entirely  eclipsed. 

6.  Brillouin  scattering  as  an  example  of  fluctuations 
about  a  hydrodynamic  nonequiljbrium  state 
Tremblay  et  al.  (1980)  considered  fluctuations  in  hydrodynamic  modes  A(k) 
governed  by  a  Langevin  matrix  equation 
9A(k, t) 

—  ~t—  +  M(k)  -A(k)  =  5(k,t)  .  (6.1) 

Notice  that  M(k)  is  the  A  operator  of  the  last  part  of  Section  A.  This  is  an 
example  of  our  transport  formalism  for  infinite  dimensional  stochastic  processes 
(y-*k)  as  considered  in  detail  in  our  1971  paper.  For  the  Langevin  forces 
correlation  we  write  as  usual 

<  ;(k,t);(k,,t’)>  =  -j-  Se(k,k')5(t- t’)  =  D(k,k’)6(t- t’)  .  (6.2) 

Solving  from  the  Langevin  equation  for  the  spectrum,  one  finds 

SA(k,k',^)  =  2[M(k)  +  i.^r1: (!$,!$')  +  2£(k,k,)[M(k)  -  lull"1  .  (6.3) 

In  a  nonequilibrium  state  F  can  be  eliminated  by  pre-  and  post-multiplying  with 
the  resolvent;  hence 

[M(k)  +  iloIjSA(k,k,,u>)[M(k)  -  i  0I  ]  =  2r(k,k')M(k) 

+  2M(k)r(k,k')  =  2D(k,k')  (6.4) 

where  we  used  the  full  noncquilibrium  second  moment  relaxation  theorem  (4.23). 
Inverting  (6.3)  we  have 

SA(k,k\„)  =  2[M(k)  +  l'nI]-1D(k,k')[|(k)  -  i,.,I] 


(6.5) 


(For  comparison  with  Tremblay  et  al.,  notice  our  S^Ck.k'.w)  is  their  2x‘J(k,k') 
with  x  =  II  xa6l;  our  Eci*  (6.5)  is  their  Eq.  (7).) 

We  now  come  to  their  specific  model.  They  consider  a  fluid  with  a 
fluctuating  stress  tensor,  in  a  small  temperature  gradient.  The  linearized 
hydrodynamic  equations  for  such  a  fluid  are 
3p(k, t) 

- —  =  -ip  c  (6.6) 

3k*v 

=  -i(k2/p)p-  [ (4 +  4q/3) /p]k? (k*v)  +  (i/p)k*S*k  (6.7) 

where  p  is  the  pressure,  p  is  the  density,  c  is  the  sound  velocity,  c,  the  bulk 
and  q  the  shear  viscosity.  There  is  no  Langevin  force  associated  with  (6.6),  but 
there  is  with  (6.7),  see  their  Eq .  (9).  Equations  (6.6)  and  (6.7),  taken  to¬ 
gether  as  a  column  matrix  equation,  are  of  the  form  (6.1),  though  the  details 

are  quite  complex.  One  can  thus  solve  for  S  (k',k",w).  They  obtain 

PP  ~  ~ 

c42T(k'  -  k")  [  2n  (k '  •  k") +  (;  -  2n/3)k' •- k""  ] 

S  (k'.^.u)  =  — - :r— ^ -  (6.3) 

PP  ~  (w:  -  c--k'-  +  i;jD  k’ •')  (  J  -  c-k"2  -  ii.D  k"') 

X  X 

where  =  (c  +  4n/3)/o.  Due  to  the  temperature  gradient  contained  in  T  =  T(r), 
the  two  Brillouin  light-scattering  peaks  obtain  an  asymmetrical  height.  One 
finds  a  difference  spectrum  6S  (u ;)  for  the  peaks  located  at  lul  =  ck,  see 

pp 

Tremblay  et  al.,  Eq.  (14).  The  experimental  values  observed  are  in  agreement 
with  the  theory,  thus  once  more  confirming  the  basic  nonequilibrium  result  (4.6) 
or  (4.23). 
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Figure  25.  Block  diagram  of  the  counting  system 
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Figure  35.  Current  noise  spectra  of  6CE6  (GE)  for  various  values 
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Figure  37.  S./S,  vs.  R  for  6AU6WC  (RCA)  at  1  KHz 


vs.  R  for  6CE6  (GE)  at  100  Hz 


Figure  39.  ST/ST  vs  R  for  6CE6  (GE)  at  1  KHz 
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Current  noise  spectra  of  6AU6VIC  (RCA)  for  two  different 
anode  voltages 
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